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1 Introduction 



In Chap. V [Magnetohydrodynamics of Protostellar Disks] it was shown that a 



sufficiently highly conducting Keplerian disk that is threaded by a weak mag- 
netic field will be subject to the magnetorotational instability (MRI) and may 
evolve into a turbulent state in which the field is strongly amplified and has 
a small-scale, disordered configuration. This turbulence has been proposed as 
O the origin of the effective viscosity that enables matter to accrete by transfer- 

-i— > ring angular momentum radially out along the disk plane. In this chapter we 

^ focus on an alternative mode of angular momentum transport that can play 

I— I an important role in protostellar disks, namely, vertical transport through the 

^ disk surfaces effected by the stresses of a large-scale, ordered magnetic field 

^ that threads the disk. 

l/-^ The possible existence of a comparatively strong, "open" magnetic field 

over much of the extent of at least some circumstellar disks around low- and 

00 intermediate-mass protostars is indicated by far-IR and submillimeter polar- 

' J ization measurements, which have discovered an ordered, hourglass-shaped 

field morphology on sub-parsec scales in several molecular clouds (e.g., Schle- 

O uning 1998; Girart ct al. 2006; Kirby 2009). The polarized radiation is at- 

tributed to thermal emission by spinning dust grains whose short axes are 
• • aligned along the magnetic field (e.g., Lazarian 2007). The detected hourglass 

_ ^ morphology arises naturally in molecular cloud cores in which a large-scale 

magnetic field provides dynamical support against the core's self-gravity (see 
;h Sect. 2.1). In this picture the field is interstellar in origin and is part of the 

^ Galactic magnetic field distribution. As discussed in Sect. 3, the inward grav- 

itational force can become dominant and the core then undergoes dynamical 
collapse to form a central protostar and a circumstellar disk. The magnetic 
field is dragged in by the infalling matter and could in principle lead to a 
large-scale "open" field configuration in the disk. 
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An ordered magnetic field that threads a disk can exert a magnetic torque 
that removes angular momentum from the interior gas. This angular mo- 
mentum can be carried away along the field lines either by torsional Alfven 
waves in a process known as magnetic braking (see Sect. 2.2) or through a 
rotating outflow in what is known as a centrifugally driven wind (CDW; see 
Sect. 4.3). These mechanisms could supplement or even entirely supplant the 
radial transport along the disk plane invoked in traditional disk models: by 
turbulent stresses as in the MRI scenario mentioned above or through grav- 
itational torques in a self-gravitating disk as described in Chap. IV [Disk 
Hydrodynamics]. In the case of radial transport, the angular momentum re- 
moved from the bulk of the matter is deposited into a small amount of gas 
at the outer edge of the disk. In the case of vertical transport, this angular 
momentum is deposited into a small fraction of the disk mass (the tenuous 
surface layers of the disk) that is removed as a CDW or else (when mag- 
netic braking operates) into the ambient medium through which the torsional 
Alfven waves propagate. The introduction of this new transport channel has 
profound implications to the structure and properties of disks in which it is 
a major contributor to the angular momentum budget and potentially also 
to the strong connection that has been found between accretion and outflow 
phenomena in young stellar objects. This is discussed in Sect. 4, where we 
also consider how to determine which of the two possible angular momen- 
tum transport modes (radial or vertical) operates at any given location in a 
magnetically threaded disk and whether these two modes can coexist. 

Of course, large-scale, ordered magnetic fields can also be produced in situ 
by a dynamo process; we consider this alternative possibility for the origin 
of the disk field in Sect. 4.2. In the case of the Sun, high-resolution observa- 
tions made at extreme ultraviolet and soft X-ray wavelengths and transformed 
into spectacular false-color images have revealed a complex web of organized 
structures that appear as loops and prominences near the stellar surface but 
simplify to a more uniform distribution further out (e.g., Balogh et al. 1995). 
There is growing evidence that Sun-like stars are already magnetically active 
in the protostellar phase and, in fact, generate fields that arc a thousand times 
stronger than that of the present-day Sun. The dynamical interaction between 
such a field and a surrounding accretion disk through which mass is being fed 
to the nascent star could have important evolutionary and observational con- 
sequences. We consider this in Sect. 5. We conclude with a summary and a 
discussion of future research directions in Sect. 6. 



2 MHD of Magnetically Threaded Disks 

Before getting into the specifics of the various topics outlined in Sect. 1, we 
present a general discussion of the dynamical properties of ordered magnetic 
fields in relation to protostellar disks and of the main methods that have been 
applied to their study. 
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2.1 Magnetic Forces 

The magnetic force per unit volume on a magnetized fluid element is given by 

where J is the current density, B is the magnetic field vector, c is the speed 
of light, and we substituted 

J = X B (2) 
47r 

(neglecting the displacement current in Ampere's Law on the assumption that 
all speeds are <C c) and used a vector identity to obtain the two terms on the 
right-hand side of Eq. (1). As was already noted in Chap. V, the first term 
represents the magnetic pressure force and the second one the force due to 
magnetic tension. To build intuition it is useful to consider the magnetic field 
lines. Because of the solenoidal (i.e., absence of monopoles) condition on the 
magnetic field, 

V • S = , (3) 

the flux ^ of "open" field lines through the disk (the integral / B ■ AS, where 
dS is a disk surface-area element) is conserved. This is a noteworthy result: it 
says that, even if the disk is resistive in its interior, the flux of magnetic field 
lines that thread it (also referred to as the poloidal flux) will not be destroyed 
(although poloidal flux could be added to, or removed from, the disk through 
its outer and inner boundaries). In the MHD picture, the field lines can be 
thought of as rubber bands that exert tension of magnitude B'^/Att directed 
along the field and pressure of magnitude B^/Sn normal to their local direction 
(e.g., Parker 1979). 

The nonmagnetic forces acting on a cloud core or a disk are the force of 
gravity (mostly self-gravity in the case of a core and central-mass gravity in 
the case of a nearly Keplerian disk), —gV(P (where g is the mass density and 
# is the gravitational potential), and the thermal pressure force —VP. There 
may be an additional force, associated with the momentum flux of MHD 
waves, in a turbulent system. Such turbulence could give rise to an effective 
viscosity, but, as is generally the case in astrophysical systems, the effects 
of microphysical shear viscosity would remain negligible (e.g., Prank et al. 
2002). Protostellar disks typically have sufficiently low temperatures that the 
dominant forces in the disk plane are gravity and magnetic tension, whereas in 
the vertical direction the thermal pressure force is always important, balancing 
the downward force of gravity or of the magnetic pressure gradient. 

These notions can be illustrated by considering the equilibrium configu- 
rations of cloud cores (e.g., Mouschovias 1976). Within the interstellar gas 
that ultimately ends up in the core the magnetic field lines arc initially nearly 
parallel to each other and the field is well coupled to the matter (with the 
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magnetic flux being, in the jargon of MHD, "frozen" into the matter). Since a 
magnetic field exerts no force parallel to itself, matter can slide relatively eas- 
ily along the field to form a flattened mass distribution. The oblate structure 
predicted by this scenario is consistent with the observed shapes of starless 
cores in the Orion giant molecular cloud (e.g., Tassis 2007). ^ It is convenient 
to describe this configuration with the help of a cylindrical coordinate system 
{r, (j), z} centered on the forming star, with the z axis aligned with the ini- 
tial magnetic field direction. A vertical hydrostatic equilibrium in which the 
(upward) thermal (and possibly turbulent) pressure force balances the (down- 
ward) gravitational force is established on a dynamical (free-fall) time scale 
and is by and large maintained throughout the subsequent evolution of the 
core, as verified by numerical simulations (e.g., Fiedler & Mouschovias 1993; 
Galli & Shu 1993), including in cases of filamentary clouds that are initially 
elongated in the field direction (e.g., Nakamura et al. 1995; Tomisaka 1996). 
Self-gravity acting in the radial direction tends to pull the field lines inward. 
The hourglass shape revealed by the polarization measurements is produced 
because these field lines remain anchored in the cloud envelope. This results 
in magnetic tension that is associated mainly with the term {Bz/4:'K)dBr/dz 
in Eq. (1) and is in complete analogy with the force exerted by a stretched 
rubber band. The field morphologies revealed by the polarization measure- 
ments are interpreted as arising from the approximate balance between this 
force and radial gravity. This interpretation is supported by H I and OH 
Zeeman measurements of the linc-of-sight field amplitude and by estimates 
of the plane-of-sky field strength using the measured dispersion in the field 
orientation (the Chandrasekhar- Fermi method), both of which typically im- 
ply roughly virialized cores (with ordered and turbulent magnetic fields con- 
tributing approximately equally to the overall support of the cloud against 
self-gravity; e.g., Ward-Thompson et al. 2007). 

In applying these ideas to the disks that form from the gravitational col- 
lapse of cores, one should bear in mind the following two points. First, the 
flux threading the disk is sufficiently strongly concentrated that the bending 
of the field lines between the midplanc (where Br = B^ = by reflection 
symmetry) and the disk surface can be large enough to make i?r,s (where the 
subscript 's' denotes the surface) comparable to B^ (which, in turn, changes 
little between the midplanc and the surface if the disk is thin) . Consequently, 
magnetic squeezing by the z-gradient of the magnetic pressure associated with 
the Br (and possibly also B^) field components can become comparable to, or 
even exceed, the downward force of gravity. This is indeed a key property of 
the wind-driving disk models discussed in Sect. 4.4. Second, the density scale 
height h in the disk typically satisfies h <^r, implying a rapid decrease of the 



^ Note, however, that many cores do not exhibit a clear oblate structure and that, 
in fact, many of the cores in Taurus are apparently prolate and may have formed 
from the fragmentation of the filamentary clouds in which they are embedded 
(e.g., Di Francesco et al. 2007). 
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density with z even as the magnetic field amplitude changes little for 2 < r. 
Therefore, magnetic forces generally dominate all other forces on scales < r 
above the disk surface, and the field there assumes a so-called "force-free" 
field configuration (J x B w 0). According to Eq. (1), in this case the mag- 
netic tension force has to balance the magnetic pressure force, which points 
outward (as increases toward the center). This implies that the field lines 
in ths region assume a vase-like, "concave in" morphology (i.e., bending to- 
ward the vertical axis), in contradistinction to the hourglass- like, "convex out" 
shape that they have inside the disk. This behavior accounts for the initial 
collimation of disk-driven MHD winds (see Sect. 4.3). 



2.2 Magnetic Braking 

Plucked rubber bands (or strings) carry waves whose phase velocity is the 
square root of the ratio of the tension to the mass density. Using again the 
analogy to magnetic field lines, one immediately obtains the phase speed of 
(transverse) Alfven waves, waz = Bz/\fiKQ (taking the background field to 
point in the z direction). When the transverse magnetic field of the wave 
points in the azimuthal direction, the corresponding (torsional) Alfven wave 
carries angular momentum. A rotating molecular cloud core or disk will twist 
the magnetic "rubber bands" that thread it. The degree of twisting fixes the 
pitch \Bff,^s/Bz\ of the field lines and is determined, in turn, by the "load" 
on the other end of the band (i.e., by the inertia of the external matter to 
which the field lines are coupled). This twisting represents a transfer of angular 
momentum from the core or disk to the external medium (subscript 'ext'). One 
can thus determine B^ ^ by equating the torque per unit area (normal to z) 
exerted on each of the two surfaces of the flattened core or disk, rB^B^^g/A-K, to 
minus the rate per unit area of angular momentum carried by the waves from 
each surface and deposited in the (initially nonrotating) ambient medium, 
—Qe^tre-xx^BVAe^ti where J7b is the angular velocity of the field line (which is 
conserved along the field under ideal-MHD conditions; see Sect. 4.3). Poloidal 
flux conservation along the field makes it possible to relate the radius r and 
field Bz in the core or disk to the corresponding quantities in the external 
medium, B^rdr = -Bz,cxt''cxtdroxt, where the ambient field (assumed to be 
uniform on large scales) can be expressed in terms of the poloidal flux ^ by 
-Bz.ext = ^/'^^Ixf then obtains (cf. Basu & Mouschovias 1994) 

Trr^ UAext 

where vb4> = rfls- The general meaning of VBtf) is discussed in Sect. 2.4; 
unless Ohm diffusivity dominates in the core/disk, it can be identified with 
the midplane angular velocity of the particles into which the field lines are 
frozen (see Eq. (65)). 
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2.3 Centrifugal Wind Driving 

As noted in Sect. 2.1, the dynamics just above the disk surface is magnet- 
ically dominated, i.e., the magnetic energy density there is larger than the 
thermal, gravitational, and kinetic energy densities of the gas. The compara- 
tively large electrical conductivity in this region implies that the poloidal (r, z) 
gas velocity is parallel to the poloidal magnetic field (see Sect. 4.3) and the 
bulk particle motions can be approximated as those of beads along rotating, 
rigid, massless wires (Henriksen & Rayburn 1971). This mechanical analogy 
is useful for deriving the criterion for the centrifugal launching of disk winds 
(Blandford & Payne 1982). We neglect thermal effects in this derivation and, 
correspondingly, regard the disk as being infinitely thin. Since the field geom- 
etry only varies on a scale ~ r in the force-free zone, very close to the disk 
surface the field lines can be regarded as being nearly straight. Considering 
thus a straight wire that intersects a Keplerian disk at a distance ro from the 
center and makes an angle 9 to the disk normal, the balance of gravitational 
and centrifugal forces along the wire implies that, in equilibrium, the effective 
potential ^cff(j/) satisfies 

d^es _ y + sinO .2 



dy (l + 2ysin6' + 2/2)3/2 



ysin2 6i-sin6l = 0, (5) 



where the dimensionless variable y = s/ro measures the distance s along the 
wire. The equilibrium is unstable when d'^'Pcs/dy'^ < (corresponding to a 
local maximum of <?eff) , which at y = occurs for 9 > 30°. Hence, centrifugal 
driving sets in if 

This wind-launching criterion plays a key role in wind-driving disk models. It is 
worth noting that, in contrast to typical stellar winds, in which the outflowing 
gas must "climb out" of basically the entire gravitational potential well at the 
stellar surface, in the case of outflows from a rotationally supported, infinitely 
thin disk the depth of the (gravitational -I- centrifugal) potential well is lower 
by a factor of 2 on account of the rotation, and gas can in principle "escape 
to infinity" without any added thermal push if there is a sufficiently strongly 
inclined, rigid channel (the magnetic field lines). Real disks, however, are not 
completely cold and therefore have a finite thickness. When this is taken into 
account, the effective potential attains a maximum at some height above the 
disk surface, and thermal pressure forces are required to lift gas from the disk 
up to that point (e.g., Ogilvie 1997). This and some other properties of CDWs 
are considered in Sect. 4.3. 



2.4 Nonideal MHD 



Unlike the approach of classical electrodynamics, in whic;li it is common to 
consider the magnetic field as being generated by current flows according to 
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the Biot-Savart law, in MHD practice it is often more illuminating to focus 
on the magnetic Held and to regard the current density as a subordinate 
quantity that is determined, through Eq. (2), by how the field is shaped by 
its interaction with matter (e.g., Parker 2007). The current, in turn, helps 
to determine the neutral-fluid-frame (denoted by a prime) electric field E' 
according to Ohm's law 

J = a ■ E' = aoE' anB X E'j_+apE'j_ , (7) 

where the conductivity has been expressed as a tensor to take account of the 
inherent anisotropy that an ordered magnetic field induces in the motions of 
charged particles. Here the subscripts || and _L denote vector components that 
are, respectively, parallel and perpendicular to the unit vector B, wherc;as 
ao, cth, and crp are, respectively, the Ohm, Hall, and Pedersen conductivity 
terms. 

Under ideal-MHD conditions, the conductivity is effectively infinite and 
the comoving electric field vanishes. This is an adequate approximation for 
describing the dynamics of a disk wind or of the medium surrounding the 
core/disk. However, within the core or disk themselves the degree of ioniza- 
tion is generally low, so finite conductivity effects must be taken into account 
in the dynamical modeling. When the conductivity is low, each charged par- 
ticle species (denoted by a subscript 'j') develops a drift velocity Wdj = Vj — v 
with respect to the neutral fluid velocity (which we approximate by the aver- 
age fluid velocity v, as appropriate for a weakly ionized medium). The drift 
velocities can be calculated from the equations of motion of these species, each 
of which is well approximated by a steady-state balance between the Lorentz 
force and the drag force F^j exerted by collisions with the neutrals, 

(£;' + ^ X b) = = mj7j^t^d,j , (8) 

where Zj is the (signed) particle charge in units of the electronic charge e, 

7j =< av >j /(wj -|- m), and < av >j is the rate coefficient for collisional 
momentum transfer between particles of mass mj and neutrals (of mass m and 
mass density g). This collisional interaction in turn exerts a force l^j„ = — -Fnj 
on the neutrals. The degree of coupling between a given charged species and 
the magnetic field is measured by the Hall parameter, defined as the ratio 
of the particle's gyrofrequency to its collision frequency Ujn = 7jQ with the 
neutrals: 

In this expression, B = sgnlB^} is the signed magnetic field amplitude, 
with the sign introduced to keep the dependence of the Hall conductivity on 
the magnetic field polarity (see Eq. (13) below). When |/3j | » 1 the coupling is 
good and the collision term can be neglected in comparison with the magnetic 
term; in this case the Lorentz force is approximately zero, corresponding to 
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the near vanishing of the electric field in the charged particle's frame. On the 
other hand, when | <C 1 the coupling is poor; in this case the magnetic term 
in Eq. (8) can be neglected. 

By summing up Eq. (8) over the particle species one obtains 

njFjn = —— . (10) 

j 

This shows expHcitly that, in a weakly ionized medium, the Lorentz force 
(which acts only on the charged particles) is transmitted to the bulk of the 
matter only through a collisional drag, which involves a relative motion be- 
tween the charged and neutral components. Therefore, if magnetic forces are 
important in such a medium, its structure is inherently not static. This is ex- 
emplified by the behavior of a magnetically supported molecular cloud core 
(e.g., Shu et al. 1987). The magnetic field that threads the core is anchored 
in the comparatively well-ionized cloud envelope, but in the core's interior 
the degree of ionization is low and the magnetic tension force is transmitted 
to the predominantly neutral gas through ion-neutral drag. Since the ions 
(taken in what follows to constitute a single species denoted by the subscript 
'i') are well coupled to the magnetic field and thus remain nearly fixed in 
space, the associated ambipolar diffusion drift entails an inward motion of 
the neutral particles toward the center of mass. If the evolution lasts longer 
than the ambipolar diffusion time (~ (-R/f A)^7i£'i) which can be inferred from 
Eqs. (14) and (18) below), the central concentration will become large enough 
to cause the core to become gravitationally unstable, and dynamical collapse 
will ensue. 

Expressing the current density in terms of the charged particles' drifts, 

•^ = Z]enjZjt;d,j (11) 
j 

(where the charged particles have number densities nj and satisfy charge neu- 
trality, X^jTijZj = 0), and using Eq. (8), one can solve for the conductivity 
tensor components in Eq. (7): 



CTH 

and 

crp 

(e.g.. Cowling 1976; Wardle & Ng 1999). Note that cth (and correspondingly 
the Hall term in Ohm's law) depends on an odd power of the magnetic field 



(12) 
(13) 

(14) 
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amplitude and can therefore have either a positive or a negative sign (reflect- 
ing the magnetic field polarity). This leads to qualitative differences in the 
behavior of the disk solutions in the Hall regime for positive and negative 
values of B^. 

The conductivity tensor formalism is useful for constructing realistic disk 
models in which the relative magnitudes of the different conductivities can 
change as a function of height even at a single radial location as a result 
of the variation in the density and in the dominant ionization mechanism 
as one moves between the disk surface and the midplane (see Sect. 4.4). It 
is nevertheless instructive to relate this formalism to the classical diffusivity 
regimes considered in Chap. V. This is best done by solving for the fluid-frame 
electric field, 

cE' = r/oV X B + 7?h(V x B) x B + tjaCV x B)j_ , (15) 
where the Ohm, Hall, and ambipolar diffusivities are given, respectively, by 

^° = 4^ ' ^''^ 
.H = ^^, (17) 

47r(T_L (Tj_ 

and 

9 

c ap 

VA = -. Vo , (18) 

47rcr_L (T_L 

with crj_ = (ct| + cr|,)^/2 (e.g., Mitchncr & Kruger 1973; Nakano et al. 2002). 
The three distinct regimes are delineated by 

(To ~ o-p > IcthI Ohm 

o-Q > IcthI > <yp Hall (19) 
(To 3> CTp 3> IcthI ambipolar . 

When ions and electrons (subscript 'e') are the only charged species, one 
has r/H = /3e Vo and ija = A (3^ r]o , with the respective Hall parameters 
given numerically by /3i = g/3e ~ 0.46(i3/ni3) (e.g., Draine et al. 1983; 
where g « 1.3 x 10^"^ Tj^^^, ni3 = nH/(10^'^ cm~^), B is measured in Gauss, 
T2 = T/(10^K), and the mean ion mass is mj — 30 mn)- In this case the 
three classical regimes correspond, respectively, to the Hall parameter ranges 
lAI < |/5e| -C 1 (Ohm), |/3i| -C 1 < |/3e| (Hall), and 1 < |A| < l^^el (ambipo- 
lar). 

Figure 1 shows the regions of ambipolar. Hall and Ohm dominance in the 
logriH ~ \ogB plane. ^ Molecular cloud cores and the outer regions (> 10 AU) 
of protostellar disks typically correspond to the ambipolar regime, as do the 



* For a complementary figure in the log riH — log T plane, see Fig. 2 of Chap. V. 
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Fig. 1. Magnetic diffusivity regimes in the lognn — logS plane for T = 280 K. 



disk surface regions at smaller radii; disk midplanes on scales ~ 1 — 10 AU 
are often dominated by Hall diffusivity, whereas Ohm diffusivity characterizes 
the midplanes of disks on scales ~ 0.1 — 1 AU. In the innermost (< 0.1 AU) 
regions of the disk, where the temperature increases above ^ 10"^ K and the 
gas becomes collisionally ionized (e.g., Gammie 1996; Li 1996a), anomalous 
Ohm diffusivity (the enhanced drag between positive and negative charge car- 
riers due to scattering off cilcictromagnctic waves generated by current-driven 
plasma instabilities) might play a role. Note that the precise extent of the dif- 
ferent regimes depends on the radial profile of the disk column density, since 
the latter determines the degree of penetration of the ionizing radiation or 
cosmic rays (see Sect. 4.4). For a given mass accretion rate, the column den- 
sity depends on the nature of the angular momentum transport mechanism. In 
particular, transport by a large-scale, ordered magnetic field is generally more 
efficient than transport by a small-scale, disordered field, ^ resulting in higher 
inflow speeds and correspondingly lower column densities and higher degrees 
of ionization in CDW-mediated accretion than in MRI-based turbulent disks. 
The behavior of the electric field is governed by Faraday's Law, 



^ This can be seen by representing the rcj) turbulent stress component as aP 
(where the constant a is typically < 1; cf. Balbus & Papaloizou 1999), so the 
ratio of the torques exerted by the ordered field and by the turbulent stress is 
{—BzB^^B/4:TT(P)){r/aH) (where () denotes vertical averaging over the disk half- 
thickness H), which is typically ^ 1. 
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(20) 



where E = E' — v x B/c is the lab-frame electric field. In view of Eq. (15), 
when the resistive term in the expression for E' can be neglected and the only 
charge carriers are ions and electrons, one can express the ambipolar and Hall 
contributions to cE as (v — Vi) x B and (vi — Vc) x B, respectively (see Konigl 
1989 and Chap. V), yielding dB/dt = V x (t>e x S), which indicates that 
in this case the field lines are frozen into the electrons (the particle species 
with the highest mobility \Z^\e/mj). It is also seen that the ideal-MHD limit 
E = —V X B/c is approached when the ion-neutral and ion-electron drift 
speeds are small in comparison with the bulk speed. It further follows that, 
in the ambipolar regime, the ions and electrons drift together relative to the 
neutrals (i.e., I^i — Ve\ <^ \vi — v\), so the field is effectively frozen also into 
the ions, whereas in the Hall regime the ions and neutrals essentially move 
together and the electrons drift relative to them {\vi — Vc\ ^ \vi — v\). More 
generally, one can define an effective velocity v b for the poloidal flux surfaces 
that applies also in the Ohm regime, when these surfaces are no longer frozen 
into any particle species (although, as noted in Sect. 2.1, they continue to 
maintain their identity). Focusing on the midplane, where B = B^, this is 
done through the relation cE — --Vb x B (cf. Umebayashi & Nakano 1986). 
As wc just observed, in the absence of Ohm diffusivity (and still considering 
only ion and electron charges) this relation is satisfied if one substitutes Vg for 
the velocity, which verifies that vb = f e in this case. However, in the Ohm 
regime the azimuthal field-line velocity 



can differ from that of the most mobile charged particle component, and 
similarly for the radial flux-surface velocity 



where the subscript '0' denotes the midplane. 

A departure from ideal MHD may lead to energy dissipation at a rate 
(per unit volume) J ■ E' . As can be inferred from Eqs. (7) and (15), both 
the Ohm and ambipolar terms in Ohm's law contribute to Joule heating, but 
not the Hall term. It may perhaps seem puzzling that energy is dissipated in 
the ambipolar regime even though the field lines remain frozen to the charged 
particles, but one can directly associate the heating in this case with the ion- 
neutral coUisional drag: J ■ E' '^it'd,i • Fin (assuming Zi 1; see Eqs. (8) 
and (11)). Joule dissipation would be the main internal heating mechanism in 
disk regions where angular momentum transport is dominated by a large-scale, 
ordered field. Ambipolar heating, in particular, could also play an important 
role in the thermal structure of disk-driven winds (see Sect. 4.3). 

Just as wc arrived at a combination of physical variables that measures the 
degree of coupling between a charged particle and the magnetic field (namely. 



(21) 



VBrfi = cE^fi/Bo , 



(22) 
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the Hall parameter; Eq. (9)), we will find it useful to identify an analogous 
measure for the neutrals. It turns out that an appropriate parameter combina- 
tion of this type for Keplerian disks, which plays a similar role in MRI-based 
systems (where the neutrals couple to a small-scale, disorderd field) and in 
wind-driving disks (where the coupling is to a large-scale, ordered field), is 
the Elsasser number 

A ^ , (23) 

/2k??_l 



where 



is the "perpendicular" magnetic difFusivity and J?k is the Keplerian angular 
velocity. This parameter is 1 and <C 1 for strong and weak neutral-field 
coupling, respectively.^ The Elsasser number is used in planetary dynamo 
theory to measure the ratio of the Lorentz to Coriolis forces. Since both the 
MRI mechanism and the field-mediated vertical angular momentum transport 
involve magnetic coupling to Keplerian rotation, it is not surprising that this 
parameter also arises in the disk context. In the three main conductivity 
regimes it reduces to 

^ = T ambipolar 

A r Hall (25) 



I. /3e A ^ Ohm , 

whore we again assumed, for simplicity, that the only charged particles are ions 
and electrons. The ambipolar limit has a clear physical meaning: it represents 
the ratio T of the Keplerian rotation time to the neutral ion momentiim- 
exchangc time. This parameter has emerged as the natural measure of the 
field matter coupling in the wind-driving disk models of Wardle & Konigl 
(1993; see Sect. 4.4) as well as in studies of the linear (e.g., Blaes & Balbus 
1994) and nonlinear (e.g., Mac Low et al. 1995; Brandenburg et al. 1995; Haw- 
ley & Stone 1998) evolution of the MRI in such disks. Indeed, since the ions 
are well coupled to the field in this limit (|/3i| ^ 1), the neutrals will be well 
coupled to the field if their momentum exchange with the charged particles 
(which is dominated by their interaction with the comparatively massive ions) 
occurs on a time scale that is short in comparison with the dynamical time 
(corresponding to T > 1). In the Hall regime A is equal to this ratio 
of time scales multiplied by This product has figured prominently in the 

® The parameter A is distinct from the Lundquist number S = vaL/tjo) and from 
the magnetic Reynolds number Rcm = VL/r]o (where V and L are characteristic 
speed and length-scale, respectively), which have been used in similar contexts in 
the literature. It was first introduced into the protostellar disk literature in this 
form (but using a different symbol) by Wardle (1999). 
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classification of wind-driving disk solutions (Konigl et al. 2010) as well as in 
linear studies of disk MRI in this regime (e.g., Wardle 1999; Balbus & Terquem 
2001). In this case, too, it has a clear physical meaning. In contradistinction 
to the ambipolar regime, the ions are not well coupled to the field in the Hall 
limit <C 1). In order for the neutrals to be well coupled to the field it 
is, therefore, not sufficient for them to be well coupled to the ions (T > 1); 
rather, the product \l3i\T must be > 1 in this case. In the Ohm regime, even 
the electrons are not well coupled to the field (|/3o| <C 1), so now the product 
^ePiT has to exceed 1 to ensure an adequate coupling of the neutrals to the 
field. The condition yl > 1 in fact characterizes both the linear (e.g., Jin 1996) 
and the nonlinear (e.g., Sano et al. 2002) behavior of the MRI in this regime, 
and the Elsasser number is again a natural parameter for classifying viable 
wind-driving disk solutions in this case (with each sub-regime corresponding 
to a distinct lower bound on A; see Konigl et al. 2010). 

2.5 Similarity Solutions 

Even the most simplified models of magnetically threaded protostellar disks 
lead to systems of nonlinear partial differential equations (PDEs). In some 
cases, such as the disk-star interaction considered in Sect. 5, it has proven nec- 
essary to carry out full-fledged numerical simulations to get definitive answers 
about the expected behavior. However, for other aspects of the problem, such 
as the disk formation (Sect. 3) and the wind-driving disk structure (Sect. 4), 
it has been possible to make analytic (or, rather, semianalytic) progress by 
looking for similarity solutions, an approach that converts the PDEs into or- 
dinary differential equations (ODEs; e.g., Landau & Lifshitz 1987). This type 
of solution does not incorporate either inner or outer radial boundaries, so 
it can be justified only at radii that are sufficiently far removed from the 
disk's actual edges. However, the basic properties revealed by these solutions 
have invariably been confirmed by simulations, and the analytic approach 
has the advantage of making the qualitative behavior more transparent and 
of making it possible to better investigate the parameter dependence of the 
solutions and to incorporate more physics and dynamic range than is yet 
possible numerically. The semianalytic and numerical approaches have thus 
played complcmcintary roles in shedding light on these questions. 

In the case of the disk-formation problem, the separate dependences on 
the spatial scale r and the time t are subsumed into a single dependence 
on the dimensionless combination r/Ct, taking advantage of the fact that 
the isothermal sound speed C is nearly uniform (on account of the efficient 
cooling by dust grains) in cloud cores and in the outer regions of circumstellar 
disks. The self-similarity of the derived solutions is expressed by the fact that 
they depend solely on the ratio r/t and not on the separate values of r and 
t. In the case of the stationary, axisymmetric wind-driving disk models, the 
separate dependences on the spatial coordinates z and r are subsumed into 
a single dependence on the dimensionless combination z/r (or, equivalently. 
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on 2;(r, ro)/?'o, where z{r,rQ) describes the shape of a magnetic field Hne that 
intersects the midplanc at Tq). In these radially self- similar models all physical 
variables scale as power laws of the spherical radius R = {z^ + r^)^/-^ and the 
differential equations depend only on the polar angle 6 = arc cot {z/r). Further 
details are given in Sees. 3 and 4, respectively. 

3 Disk Formation and Early Evolution 

Rotationally supported circumstellar disks evidently originate in the collapse 
of self-gravitating, rotating, molecular cloud cores. Molecular line observa- 
tions (e.g., Goodman et al. 1993; Kane & Clemens 1997) have established 
that a majority of dense (> 10* cm~^) cloud cores show evidence of rota- 
tion, with angular velocities ~ 3xl0~^^ — 10"^'^ s~^ that tend to be uni- 
form on scales of ~ 0.1 pc, and with specific angular momenta in the range 
~ 4x10^° — 3x10^^ cm^ s~^. The cores can transfer angular momentum to the 
ambient gas by magnetic braking (Sect. 2.2), and this mechanism also acts 
to align their angular momentum vectors with the local large-scale magnetic 
field (e.g., Machida et al. 2006). This alignment occurs on a dynamical time 
scale and hence can be achieved even in cores whose lifetimes are not much 
longer than that (as in certain models of the turbulent ISM; e.g., Elmegreen 
2000). The dynamical collapse might occur as a result of mass rearrangement 
in the core on the ambipolar diffusion time (e.g., Mouschovias et al. 2006; sec 
Sect. 2.4) or sooner if the core is close to the critical mass for collapse (the 
effective Jeans mass) from the start (e.g., Elmegreen 2007). In this section we 
consider the core collapse problem in the context of angular momentum trans- 
port by a large-scale, ordered magnetic field; an alternative scenario involving 
gravitational torques is discussed in Chap. IV. 

Once dynamical collapse is initiated and a core goes into a near-free-fall 
state, the specific angular momentum is expected to be approximately con- 
served, resulting in a progressive increase in the centrifugal force that eventu- 
ally halts the collapse and gives rise to a rotationally supported disk on scales 
~ 10^ AU. These expectations are consistent with the results of molecular- 
line interferometric observations of contracting cloud cores (e.g., Ohashi et 
al. 1997; Belloche et al. 2002). In this picture, the disk rotation axis should 
be aligned with the direction of the large-scale magnetic field that threads 
the cloud. Observations have not yet yielded a clear-cut answer to whether 
this is indeed the case in reality (e.g., Menard & Duchenc 2004; Vink et al. 
2005), and it is conceivable that the field in some cases is too weak to align 
the core's rotation axis or even control its contraction, or that additional pro- 
cc!ss(^s (such as fragmentation, disk warping, etc.) play a role. In what follows 
we nevertheless continue to pursue the implications of the basic magnetically 
supported cloud picture. 
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3.1 Modeling Framework 

Since the gravitational collapse time is much shorter than the local ambipolar 
diffusion time in the core, the magnetic field lines at first move in with the 
infalling matter. However, once the central mass begins to grow, ambipolar 
diffusion becomes important within the gravitational "sphere of influence" of 
the central mass (Ciolek &: Konigl 1998; Contopoulos et al. 1998). When the 
incoming matter enters this region, it decouples from the field and continues 
moving inward. The decoupling front, in turn, moves outward and steepens 
into an ambipolar diffusion shock. The existence of this C-type MHD shock 
was first predicted by Li & McKee (1996). The transition from a nearly freely 
falling, collapsing core to a quasi-stationary, rotationally supported disk in- 
volves a strong deceleration in a centrifugal shock. This shock typically occurs 
at a smaller radius than the ambipolar-diffusion shock and is hydrodynamic, 
rather than hydromagnetic, in nature. 

Based on the arguments presented in Sect. 2.1 it should be a good approx- 
imation to assume that the gas rapidly establishes force equilibrium along the 
field and therefore to consider only motions in the radial direction.'^ One can 
obtain semianalytic solutions for this effectively ID time-dependent problem 
by postulating r — t self-similarity, with a similarity variable 



(see Sect. 2.5). This modeling approach is motivated by the fact that core 
collapse is a multiscale problem, which is expected to assume a self-similar 
form away from the outer and inner boundaries and not too close to the onset 
time (e.g., Penston 1969; Larson 1969; Hunter 1977; Shu 1977). This has been 
verified by numerical and semianalytic treatments of restricted core-collapse 
problems - with/without rotation and with/without magnetic fields. Although 
the constancy of the isothermal soimd speed C that underlies the ansatz (26) 
is not strictly maintained throughout the solution domain (in particular, C 
scales roughly as r~^/^ in large portions of the disk that forms around the 
central star) . this is of little consequence to the results since thermal stresses 
do not play a major role in the dynamics of the collapsing core. For a typical 
sound speed C = 0.19 kms'S x = l^ {400, 4000} AU at t = {10^ lO^jyr. 

Krasnopolsky & Konigl (2002, hereafter KK02) constructed self-similar 
solutions of rotating magnetic molecular cloud cores that are subject to am- 
bipolar diffusion. These solutions reveal many of the basic features of star and 
disk formation in the core-collapse scenario and arc; discussed in the remain- 
der of this section. To incorporate ambipolar diffusion into the self-similarity 

^ In reality, mass can also be added to the system from the polar directions. Numer- 
ical simulations of axisymmetric collapse in which new mass is added only through 
vertical infall have tended to produce disk-to-star mass ratios ~ 1, much higher 
than typically observed. A better agreement with observations may, however, be 
obtained when nona:xisymmetric density perturbations and resultant gravitational 
torques are included in the calculations (e.g., Vorobyov & Basu 2007). 



r 



(26) 
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formulation it is necessary to assume that the ion density scales as the square 
root of the neutral density: Qi = ICg^^"^. As discussed in KK02, this should be 
a good approximation for the core-collapse problem: it applies on both ends of 
a density range spanning ~ 8 orders of magnitude, which corresponds roughly 
to radial scales ~ 10 — 10* AU, with /C varying by only ^ 1 order of magnitude 
across this interval. 

To allow mass to accumulate at the center in a ID rotational collapse, 
an angular momentum transport mechanism must be present. KK02 assumed 
that vertical transport through magnetic braking continues to operate also 
during the collapse phase of the core evolution. To incorporate this mecha- 
nism into the self-similar model it is necessary to assume that WAcxt , the Alfvcn 
speed in the external medium, is a constant. * KK02 verified that, in their 
derived solutions, magnetic braking indeed dominates the most likely alter- 
native angiilar-momentum transport mechanisms — MRI-induccd turbulence 
and gravitational torques. However, they also found that angular momentum 
transport by a CDW arises naturally (and may dominate) in the Keplerian 
disk that forms in their fiducial solution (see Sect. 3.4). 



3.2 Basic Equations 

The mass and momentum conservation relations in their differential form are 

given by: 

Mass 

do 1 d , , 9 , , 

Radial Momentum 



at or Or r 47r oz 

a fBi\ 1 a , „ ,0 du 



Angular Momentum 



dr V Stt / 87rr^ dr dz 



{rB^f-gv,^, (28) 



g d gvr d B^ dB,i, Br d d 

rdt + —d-r = 4;^^ + 4^^^^;^ (^^^^ - ' (29) 

Vertical Hydrostatic Equilibrium 

^'1 = ^^5.-^1^ + ^!+^^, (30) 




A nearly constant value UAext ~ 1 km s is, in fact, indicated in molecular clouds 
in the density range ~ 10^ - 10^ cm"^ (e.g., Crutcher 1999). 



Chapter VI: The Effects of Large-Scale Magnetic Fields 



17 



where gr and gz are, respectively, the radial and vertical components of the 
gravitational acceleration. We have not yet dropped the terms that involve 
the vertical velocity component (except in the vertical momentum equation, 
which we assume takes on its hydrostatic form). 

We now integrate these equations over the core/disk thickness 2H. We use 
the thin-disk approximation {H{r) <C r) and assume that the density, radial 
velocity, azimuthal velocity, and radial gravity are constant with height, and 
that so is also (except when it is explicitly differentiated with respect to 
z, in which case we substitute dB^/dz = —r~^{d/dr){rBr) from V ■ B = 
0). We also assume Br{r,z) = Br^sir) [z/H{r)] (and similarly for B^), and, 
after deriving expressions that are valid to order (H/r)"^, further simplify by 
dropping all terms 0{H/r) except in the combination [Bj-s — H{dBz/dr)\ 
(see KK02 for details). We thus obtain the vertically integrated conservation 
equations: 
Mass 

^ + -^(^^M = -2^^, (31) 

Radial Momentum 

dvr ^ dvr C^dS Bz ( rrdBA 

-dt^''^^=''-U^^2^ i^- - J + r3 ' (32) 

Angular Momentum 



dJ dJ rBzBd, 



dt dr 27rr ' 
Vertical Hydrostatic Equilibrium 



(33) 



W - 2 + 2r3 + 8^^ + ^'•■^ - ^'•■^ j ' ^^^^ 

where S = 2gH is the surface mass density, J = ru^ is the specific angular 
momentum of the matter, G is the gravitational constant, and is the mass 
of the central protostar. In the integrated equations we have implemented 
the ID flow approximation by setting ~ 0, but we retained the term on 
the right-hand side of Eq. (31) to allow for mass loss through a disk wind 
(at a rate Mw). Such a wind could carry angular momentum, but we did not 
include vertical particle angular momentum transport in Eq. (33) since most 
of a disk wind's angular momentum initially resides in the magnetic field (see 
Sect. 4.3). In any case, we proceed to solve the equations by assuming at first 
that no wind is present and that the only mechanism of angular momentum 
transport is magnetic braking. 

The dominant charge carriers in the pre-coUapse core are ions and elec- 
trons. Adopting this composition, we approximate the ion equation of motion 
in the ambipolar-diffusion limit by 



18 



Arieh Konigl and Raquel Salmeron 



gi^niVd^^{yy-B)xB, (35) 

where v^i is the neutral-ion collision frequency (see Eqs. (8) and (10)). This 
relation yields the components of the drift velocity: 

Magnetic braking is incorporated via Eq. (4), in which we identify 
with Vi^cj, (see Sect. 2.4). Expressing Vi^cf, in terms of ii^ and Vd,cj, (Eq. (36)) and 
imposing a cap {6 < 1) on \B^^s/Bz\ (to account for the possible development 
of a kink instability above the disk surface), one gets 



<!' V4, f ^Bz 



SB, 



(38) 



The flux conservation relation d\^/dt = —2'KrvBrBz (obtained from 
Eqs. (20) and (22)), which describes the advection of poloidal flux by the 
infalling matter, can be written in this limit as 

— = = -2-Kr {Vr + Ud,r) , (39) 

with Wd,r given by Eq. (37). The surface value of B^. that appears in the latter 

equation can be determined in the limit of a potential field (V x B = 0) 
outside an infinitely thin disk from an r integral of the midplane value of B, '■ 



/■OO /"OO 

Br,s= dkJiikr) drV[B^(r') - Bref] Jo(fcr') 
Jo Jo 



(40) 



(e.g., Ciolek & Mouschovias 1993), where Jq and Ji are Bessel functions of 
order and 1, respectively, and Brcf is the uniform ambient field at "infin- 
ity" (which is henceforth neglected). Although the current-free limit of the 
force- free medium outside the core/disk is not exact, as the presence of a 
component, which typically involves a distributed poloidal electric current, is 
required for angular momentum transport above and below the core/disk sur- 
faces, the magnitude of \B^^s/Bz\ (as deduced from Eq. (38)) is typically small 
at the start of the collapse (e.g., Basu & Mouschovias 1994). Furthermore, even 
after this ratio increases to ~ 1 as matter and field become concentrated near 
the center, the magnitude of the enclosed (poloidal) current (7 = {c/2)r\B^\) 
increases relatively slowly with decreasing r (scaling as r^^/^ in the circum- 
stellar disk; see Eq. (55) below), so overall we expect the expression (40) to 
be a fair approximation. A similar integral (with S{r) replacing B,) can be 
written for gr (Toomre 1963). These integrals can be approximated by their 
monopole terms 
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(cf. Li & Shu 1997), where M{r,t) is the mass enclosed within a radius r at 
time t. When Bz{r) (or S{r)) scales as r"'', one still obtains the nionopole 
expression for Br,s (or Qr), but with a g-dependent coefScient 0(1) (Ciolek & 
Konigl 1998). 



3.3 Self-Similarity Formulation 

The various physical quantities of the problem can be expressed as dimension- 
less functions of the similarity variable x (Eq. (26)) in the following fashion: 

H{r, t) = Ct h{x) , E{r, t) = {C/2nGt) a{x) , (43) 

Vr{r,t) = Cu{x) , v^{r,t) = Cw{x) , (44) 

gr{r,t) = {C/t)g{x), J{r,t) = CH 3{x) , (45) 

M(r, t) = {CH/G) m{x) , M^{r, t) = (C^/G) m{x) , (46) 

B{r, t) = {C/G^'H) b{x) , ^^(r, t) = {2TrCh/G^^^)iP{x) (47) 

where is the mass accretion rate. 

The vertical force-balance relation (34) yields a quadratic equation for the 
normalized disk half-thickness h, 

- ^r/^) + + hi, + a^)h-2a = Q, (48) 

whose solution is 

where m, = m* — x^br,{dbz/dx)/a and a = a + {b^, + b^ J/a. 

The initial (t = 0) conditions, which in the self-similar model also represent 
the o'ider asymptotic (r — >• oo) values, correspond to a collapsing core just 
before it forms a central point mass. Based on previous analytic and numerical 
work, KK02 adopted 

A a 

<j ^ —, bz ^ , u — >■ Uoo, w Woo as a; — >■ oo , (50) 

X l^oo 

with ^ = 3, /Zoo = 2.9, = —1- Prom the constituent equations one can 
also derive the inner asymptotic behavior (corresponding to r ^ at a fixed 
t) in the Keplerian disk: 



h = 



ax 
2m» 
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m = m 



—u 



1/2 1/2 

(m*/cri)a;^/^ , 

(2r?/3J)(2m,)V^ -3/2 
[1 + (2t/3J)-2]i/2'^ 

(4/3)6, , 

{2/[l + (2T/3(5)2]m4V2a.3/2 ^ 



(51) 
(52) 
(53) 



bz = -b4,,s/S 
br,s = tp/x^ 
h 



(54) 
(55) 
(56) 
(57) 



where r = {iirGgy/^ /u^i. 

The equations arc solved as a boundary vahic problem using the above 
asymptotic relations. A solution is determined by the values of the four model 
parameters: t, 6 = \B^s/Bz\, Woo = v^{t = 0)/C, and a = C/^Aext- The 
scaling parameter to* for the central mass and the mass accretion rate is 
obtained as an eigenvalue of the problem. 

3.4 Self-Similar Collapse Solutions 

Figure 2 shows a fiducial solution corresponding to the parameter combination 

T = 1, 5 = 1, w(0) = 0.73, and a = 0.8 (which yield to* = 4.7). In this case 
the initial rotation is not very fast and the braking is moderate, leading to 
the formation of a disk (with outer boundary at the centrifugal shock radius 
Xc = 1.3 X 10~2) within the ambipolar-diffusion (AD) region (delimited by 
the AD shock radius Xa = 0.41 « 30 Xc)- One can distinguish the following 
main flow regimes: 

• Outer region {x > x^)- ideal-MHD infall. 

• AD shock: resolved as a continuous transition (but may in some cases 
contain a viscous subshock); KK02 estimated Xi, « \f2T j \x^. 

• Ambipolar diffusion-dominated infall (sc < a; < Xs): near free fall con- 
trolled by the central star's gravity. 

• Centrifugal shock: its location depends sensitively on the diffusivity param- 
eter T, which affects the amount of magnetic braking for x < a;a; KK02 
estimated x^ « (m»w^/A2)cxp{-(2^/2m,//ioo)"^^^'''^^/^}. 

• Keplerian disk (x < Xc): asymptotic behavior (Eqs. (51) (57)) is ap- 
proached after a transition zone representing a massive ring (of width 
~ 0.1 Xc and mass ~ 8% of the disk mass within Xc, which itself is < 5% 

of TO*). 

The inner asymptotic solution implies that, at any given time, the disk 
satisfies MJ^r) = const, Eir) cx r'^/^, B oc r'^/^, and = 4/3.^ If 

^ The surface density scaling is the same as that inferred for the "minimum mass" 
solar nebula (e.g., Weidenschilling 1977). Note, however, that this scaling is also 
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Fig. 2. Behavior of key normalized variables in the fiducial self-similar disk formar- 
tion solution as a function of the similarity variable x (see Eqs. (43)-(46)). 



the derived power-law dependence of Bz on r is used in Eq. (40) one obtains 
Br,s/Bz = 1.428, slightly less than the result found with the monopole ap- 
proximation (41), but still representing a strong bending of the field lines (by 
~ 55° from the normal) that is significantly larger than the minimum of 30° 
required for launching a centrifugally driven wind from a "cold" Keplerian 
disk (Eq. (6)). This implies that protostcUar disks formed in this fashion are 
likely to drive disk outflows over much of their radial extents, at least during 
early times when they still accumulate mass from the collapsing core. Inter- 
estingly, the steady- state, radially self-similar CDW solution of Blandford & 
Payne (1982) also yields a radial magnetic fleld scaling oc r~^/'^. As discussed 
by KK02, this suggests that angular momentum transport by a CDW can be 
formally incorporated into the disk-formation solution, modifying some of the 
details (such as the mass accretion rate onto the star, which could be reduced 
by up to a factor of ~ 3) although not its basic properties. When interpreted 
as a wind-driving disk model, the asymptotic solution corresponds to a weakly 
coupled disk configuration (using the nomenclature introduced in Sect. 4.4). 

Despite being based on a number of simplifying assumptions, the fiducial 
solution demonstrates that vertical angular momentum transport along inter- 



predicted for a self-similar Keplerian disk with an "a viscosity" (Tsuribe 1999) 
as well as in certain models in which gravitational torques dominate the angular 
momentum transport (e.g., Lin & Pringle 1987; Voroboyov & Basu 2007). 



22 



Arieh Konigl and Raquel Salmeron 



stellar magnetic field lines can be sufficiently efficient to allow most of the 
mass of a collapsing molecular cloud core to end up (with effectively no angu- 
lar momentum) at the center, with the central mass dominating the dynamics 
well beyond the outer edge of the disk even as the inflow is still in progress. 
The vertical transport is thiis seen to resolve the so-called angular momen- 
tum problem in star formation (although the exact value of the protostar's 
angular momentum is determined by processes near the stellar surface that 
are not included in this model; see Sect. 5). To the extent that self-similarity 
is a good approximation to the situation in such cores, it is conceivable that T 
Tauri (Class II) protostellar systems, whose disk masses are typically inferred 
to be < 10% of the central mass, have had a similarly low disk-to-star mass 
ratio also during their earlier (Class and Class I) evolutionary phases. This 
possibility remains to be tested by observations. 

The fiducial solution also reveals that the AD shock, even though it is 
located well outside the region where the centrifugal force becomes impor- 
tant, helps to enhance the efficiency of angular momentum transport through 
the magnetic field amplification that it induces. The revitalization of ambipo- 
lar diffusion behind the AD shock in turn goes a long way toward resolving 
the magnetic flux problem in star formation (the several-orders-of-magnitude 
discrepancy between the empirical upper limit on the magnetic fiux of a pro- 
tostar and the flux associated with the corresponding mass in the pre-coUapse 
core), although it is conceivable that Ohm diffusivity in the innermost regions 
of the disk also plays an important role in this process (e.g., Shu et al. 2006; 
Tassis & Mouschovias 2007). The details of the flux detachment outside the 
Ohm regime can be modifled if one takes account of the fact that the flux 
is strictly frozen into the electrons but not necessarily into the ions (i.e., if 
one includes also the Hall term in Ohm's Law; see Sect. 2.4). In particular, 
Tassis & Mouschovias (2005) found that a quasi-periodic series of outward- 
propagating shocks may develop in this case. 

The shocks that form in the collapsing core, and in particular the cen- 
trifugal shock (which is the strongest), will process the incoming material and 
may have implications to the composition of protoplanctary disks (e.g., the 
annealing of silicate dust; see Marker & Desch 2002 and Chap. VIII [Dust 
Processing and Mineralogy in Protoplanetary Accretion Disks]). Note in this 
connection that the shock velocity Vgh relative to the intercepted matter has 
roughly the free-fall magnitude oc {M/ry^^, which is constant for any given 
value (including, in particular, Xc) of the similarity variable. The postshock 

1 /2 

temperature (which scales as v^^i ) '^iU thus not vary with time if the inflow 
remains self-similar; the postshock density, on the other hand, will decrease 
with time as the shock moves to larger radii. 

By modifying the model parameters, one can study the range of possible 
behaviors in real systems. Figure 3 shows two limiting cases, which bracket 
the fiducial solution (cf. Fig. 2). 
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Fig. 3. Behavior of key normalized variables in representative fast-rotation (left) 
and strong-braking (right) solutions of the collapse of a rotating, magnetic, molec- 
ular cloud core as a function of the similarity variable x (see Eqs. (43)-(47)). 



The fast rotation case differs from the fiducial sohition primarily in having a 
large initial-rotation parameter (u'(O) = 1.5). It has the following distinguish- 
ing features: 

• The centrifugal shock is located within the self-gravity-dominated (and 
idcal-MHD) region; a back-flowing region is present just behind the shock. 

• The central mass is comparatively small (m* = 0.5), giving rise to a non- 
Keplerian outer disk region. 

• The ideal-MHD/ambipolar-diffusion transition occurs behind the centrifu- 
gal shock and is gradual rather than sharp. 

The strong braking case (r = 0.5, 6 = 10, w(0) = 1, a = 10, yielding m* = 
5.9) is characterized by large values of the braking parameters d and a. It is 
distinguished by having 

• no centrifugal shock (or circumstellar disk); essentially all the angular 
momentum is removed well before the inflowing gas reaches the center, and 
the X ^ behavior resembles that of the nonrotating collapse solution of 
Contopoulos et al. (1998). 

This case may be relevant to the interpretation of slowly rotating young stars 
that show no evidence of a circumstellar disk (e.g., Stassun et al. 1999, 2001), 
which would seem puzzling if all protostars were born with significant rotation 
that could only be reduced as a result of their interaction with a disk (see 
Sect. 5). In the strong-braking interpretation, both the slow rotation and 
the absence of a disk have the same cause. It is also worth noting in this 
connection that the physical conditions that underlie the fiducial solution 
presented above — namely, a highly flattened, magnetized core whose braking 
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is controlled by the inertia of a comparatively low-density ambient gas — may 
not be generally applicable. In particular, it has been argued (see Mellon & 
Li 2009 and references therein) that magnetized molecular clouds may retain 
a significant amount of mass outside the collapsing core that would dominate 
the braking process and could place the core in the strong-braking regime 
where disk formation would be suppressed. Although there may be ways to 
alleviate or even circumvent this problem (e.g., Hennebelle & Ciardi 2009; 
Dufhn & Pudritz 2009) , there is obviously a need for additional observational 
and theoretical work to clarify the physical conditions that characterize core 
collapse in real clouds. 

While the scmianalytic self-similarity model described in this section is 
useful for identifying important features of disk formation from the collapse of 
rotating cores, the intricacies of this process can only be studied numerically. 
Newly developed 3D, nonideal-MHD, nested-grid codes (e.g., Machida et al. 
2007) are already starting to provide new insights into this problem. 

4 Wind-Driving Protostellar Disks 

In this section we discuss wind-driving disk models. We begin by highlighting 
the observational evidence for a close link between outflows and accretion disks 
in protostellar systems and the interpretation of this connection in terms of a 
large-scale magnetic field that mediates the accretion and outflow processes. 
We next outline the theory of centrifugally driven winds, describe the incor- 
poration of CDWs into equilibrium disk models, and conclude by addressing 
the stability of such disk/wind configurations. 

4.1 The Disk— Wind Connection 

We give a brief survey of the observational findings on protostellar outflows 
and their connection to accretion disks. Besides the explicit references cited, 
the reader may consult some of the many review articles written on this topic 
over the years (e.g., Cabrit 2007; Calvet et al. 2000; Dutrey et al. 2007; Konigl 
& Pudritz 2000; MiUan-Gabet et al. 2007; Najita et al. 2007; Ray et al. 2007; 
Watson et al. 2007). 

Bipolar Outflows and Jets 

Bipolar molecular outflows and narrow atomic (but sometimes also molecu- 
lar) jets are ubiquitous phenomena in protostars, with nearly 1000 collimated 
outflows of all sorts already known (e.g.. Bally et al. 2007). The bipolar lobes 
are usually understood to represent ambient molecular material that has been 
swept up by the much faster, highly supersonic jets and winds that emanate 
from the central star/disk system (e.g., Arce et al. 2007). Jets associated with 
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low-bolometric-luminosity (I/boi < 10^ Lq) protostars have velocities in the 
range ^ 150 — 400 km s^^, large (> 20) Mach numbers, and inferred mass 
outflow rates ~ 10~^ — 10"'' M0yr~^. They are coUimated on scales of a 
few 10s of AU and exhibit opening angles as small as ~ 3 — 5° on scales of 
10^ — 10** AU. High-resolution observations of optically visible jets from clas- 
sical T Tauri stars (CTTSs) reveal an onion-like morphology, with the regions 
closer to the axis having higher velocities and excitations and appearing to 
be more collimatcd. Detailed optical/ncar-IR spectral diagnostic techniques 
have been developed and applied to classical (Herbig-Haro) protostellar jets, 
making it possible to directly estimate the neutral densities in the forbidden- 
line emission regions. Recent results (Podio ct al. 2006) yield an average mass 
outflow rate of 5 x lO~*M0yr~^, markedly higher than previous estimates. 
Outflows have also been detected by optical observations of intermediate-mass 
(^2— 10 M0) Herbig Ac/Be stars and other high-luminosity sources. The jet 
speeds and mass outflow rates in these protostars are larger by a factor ~ 2 — 3 
and ~ 10 — 100, respectively, than in the low-iboi objects. 

The Link to Accretion Disks 

Strong correlations are found between the presence of outflow signatures (P- 
Cygni line proflles, forbidden line emission, thermal radio radiation, well- 
developed molecular lobes) and accretion diagnostics (UV, IR, and millimeter 
emission excesses, inverse P-Cygni line profiles) in T Tami stars (e.g., Har- 
tigan et al. 1995). Such correlations evidently extend smoothly to protostars 
with masses of ~ 10 Mq. A related finding is that the apparent decline in 
outflow activity with stellar age follows a similar trend exhibited by disk fre- 
quency and inferred mass accretion rate. In addition, correlations of the type 
M oc L^^^j (with q ~ 0.6 — 0.7) have been found in both low-Lboi and high-iboi 
protostars for mass accretion rates and for mass outflow rates in ionized jets 
as well as in bipolar molecular lobes. Furthermore, CTTS-like accretion and 
outflow phenomena have now been detected also in very-low-mass stars and 
brown dwarfs (e.g., Mohanty et al. 2005). 

These flndings strongly suggest that outflows are powered by accretion and 
that the same basic physical mechanism operates in both low- (down to nearly 
the planetary mass limit) and intc!rnicdiate-mass protostars, and possibly also 
in some higher-mass objects (e.g.. Shepherd 2005). The accretion proceeds 
through circumstellar disks, which can be directly probed by means of high 
spectral- and spatial-resolution (in particular, interferometric) observations 
at submillimeter, millimeter, mid-IR, near-IR, and optical wavelengths. The 
disks appear to be rotationally supported (for r < 100 AU), and when the 
rotation law can be determined, it is usually consistent with being Keplerian 
{V(f, oc r~^/^). Recent data indicate that M^isk ^ at least up to ~ 
20 M©. 

Strong evidence for a disk origin for the observed outflows is available 
for FU Orionis outbursts in rapidly accreting young protostars (e.g.. Hart- 
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mann & Kenyon 1996). It is inferred that the emission during an outburst (of 
typical duration ^ 10^ yr) originates in a rotating disk and that the outflow 
represents a wind that accelerates from the disk surface (with M^jM^ ~ 0.1; 
Ma ~ 10~* Mq yr~^). It has been suggested (Hartmann 1997) that most of the 
mass accumulation and ejection in low-mass protostars occurs during recur- 
rent outbursts of this type. Interestingly, estimates in CTTSs (e.g., Kurosawa 
et al. 2006; Ray et al. 2007) indicate that M^/Ma has a similar value (i.e., 
0.1) also during the quiescent phases of these protostars. 

4.2 Magnetic Driving of Protostellar Outflow^s 

The most widely accepted explanation of protostellar outflows is that they 
tap the rotational kinetic energy of the disk (and/or central object) and are 
accelerated centrifugally from the disk (or stellar) surface by the stress of a 
large-scale, ordered magnetic field that threads the source. In this subsection 
we review the arguments that have led to this picture and list a few alterna- 
tives for the origin of the field. 

Outflow Driving Forces 

The momentum discharges inferred from observations of protostellar jets are 
compatible with the values deduced for the bipolar molecular outflows but are 
typically a factor ~ 10^ — 10^ higher than the radiation-pressure thrust iboi/c 
prodiiced by the stellar luminosity, ruling out radiative acceleration as a dom- 
inant driving mechanism in low/intermediate-mass protostars (e.g., Richer et 
al. 2000). Radiative eff'ects could however be important in driving photoevap- 
orative disk outflows, particiflarly in high-Lboi systems (e.g., HoUenbach et 
al. 1994; see Chap. VII [The Dispersal of Disks around Young Stars]). Disk 
heating by a luminous protostar (such as a Herbig Be star) could potentially 
also induce a line-driven wind from the inner disk (Drew et al. 1998). An- 
other effect is radiation pressure on the dusty outer regions of a disk wind, 
which could act to decollimate the streamlines even in comparatively low-iboi 
systems (see Sect. 4.3 and cf. Konigl & Kartje 1994). 

Thermal pressure acceleration is commonly discounted as a dominant 
mechanism since the requisite high (effectively virial) temperatures are not 
generally observed at the base of the flow. It has, however, been suggested 
that, under suitable conditions, the thermal energy released in shocks at the 
boundary layer between the disk and the star could be efficiently converted 
into outflow kinetic energy (Torbett 1984; Soker & Regev 2003). Even if this 
is not a dominant mechanism, thermal pressure effects are nevertheless im- 
portant in the mass loading of hydromagnetic winds (see Sect. 4.3) and could 
potentially also play a significant role in the initial acceleration of disk out- 
flows of this type (e.g., Pesenti et al. 2004). 

A general result for protostellar jets can be obtained by combining 
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1. Mjetl-jet ~ 10^ - 10^ ibol/c 

3. ujot - 10-3 c 

4. Lbol - iacc GM.MJR,) 

(where i?* is the stellar radius), which implies that, on average, the jet kinetic 
luminosity is related to the released accretion luminosity by Ljet ~ 0.1 — 
1 L^cc-^^ Such a high ejection efficiency is most naturally understood if the jets 
are driven hydromagnetically (see Eqs. (68) and (70) below). This mechanism 
also provides a natural explanation of the strong coUimation exhibited by 
protostellar jets (see Sect. 4.3). 

Origin of Large-Scale Disk Magnetic Field 

Perhaps the most likely origin of a large-scale, open magnetic field that can 
launch a CDW from a protostellar disk is the interstellar field that supports 
the natal molecular cloud core. We have already considered the observational 
evidence for such a core-threading field (Sect. 1), its dynamical effect on the 
core (Sect. 2.1), and how it is advected inward by the inflowing matter during 
gravitational collapse (Sect. 3). However, to complete the argument it is still 
necessary to demonstrate that, once a rotationally supported disk is estab- 
lished, magnetic flux that is brought in to its outer edge remains sufficiently 
well coupled to the inflowing matter to be distributed along the disk plane. 
A potential problem arises, however, in turbulent disks in which only radial 
angular momentum transport is present. In particular, in a disk of charac- 
teristic radiiis R and half-thickness _ff <C i? in which angular momentum 
transport is due to an effective turbulent viscosity I'turb) inward dragging of 
magnetic field lines is possible only if the effective turbulent magnetic diffu- 
sivity r^turb satisfies rj < (H/R) v,^^ which may not occur naturally in such a 
system (Lubow et al. 1994a). The introduction of vertical angular momen- 
tum transport along a large-scale magnetic field offers a straightforward way 
out of this potential dilemma by decoupling these two processes: magnetic 
diffusivity is determined by the ionization state and composition of the disk 



Note that in protostars that are beyond the main accretion phase Lboi can be ex- 
pected to exceed Lace (e.g., Tilling et al. 2008), which would increase the inferred 
value of Ljet/I/acc in these sources. 

This follows from a comparison of the mass inflow speed, ~ vtuib/R, as inferred 
from the angular momentum conservation equation, with the field diffusion speed 
relative to the gas, ~ qturh/H, as inferred from the flux conservation equation. 
Recent numerical simulations have indicated that the turbulent rj is, in fact, 
of the order of the turbulent v in disks where the turbulence derives from the 
magnetorotational instability (Guan & Gammie 2009; Lesur & Longaretti 2009). 
Note, however, that efficient inward advection of magnetic flux may nevertheless 
be possible in this case on account of the expected suppression of the instability 
near the disk surfaces (e.g., Rothstein & Lovelace 2008). 
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(or else by the local turbulence, if one is present) and no longer has the same 
underlying physical mechanism as the angular momentum transport. The disk 
formation models presented in Sect. 3 provide concrete examples of how open 
magnetic field lines can be self-consistently incorporated into a disk in which 
they constitute the dominant angular momentum transport channel (see also 
Spruit & Uzdensky 2005). 

An alternative possibility is that the outflow is driven by the stellar 
dynamo-generated field. This could happen cither along field lines that have 
been effectively severed after they penetrated the disk, as in the X-wind sce- 
nario (e.g., Shu et al. 2000), or along field lines that are still attached to the 
star. However, even in the latter case the outflow is envisioned to result from 
an interaction between the stellar magnetic field and the disk (see Sect. 5). It 
has also been proposed that a large-scale, open field could be generated by a 
disk dynamo (e.g.. Tout & Pringle 1996; von Rekowski et al. 2003; Blackman & 
Tan 2004; Pudritz et al. 2007; Uzdensky & Goodman 2008). In this connection 
it is worth noting that, even if the disk dynamo generates small-scale, closed 
magnetic loops, it is conceivable, extrapolating from the situation on the Sun, 
that some of these loops could be dynamically extended to sufficiently long 
distances (in particular, beyond the respective Alfven surfaces) for them to 
become effectively open (Blandford & Payne 1982). 

Direct evidence for the presence of magnetic flelds in either the outflow 
or the disk has been scant, but this may not be surprising in view of the 
considerable observational challenges. Among the few actual measurements 
are the strong circular polarization that was detected on scales of ~ 20 AU in 
T Tau S (Ray et al. 1997), which was interpreted as a field of at least several 
gauss that was advected from the origin by the associated outflow, '^'^ and 
the meteoritic evidence that points to the presence in the protosolar nebula 
of a ~ IG field at r ~ 3 AU (Levy & Sonett 1978). In view of the strong 
indications for disk-launched outflows in FU Orionis systems, it is also worth 
mentioning the Zeeman-signature least-square deconvolution measurement in 
FU Ori (Donati et al. 2005), which was interpreted as a ~ IkG field (with 
l-B^I ^ Bz/2) originating on scales of ~ 0.05 AU in the associated disk (with 
the direction of being consistent with its origin in shearing of the vertical 
field by the disk differential rotation). Further support for a magnetic disk- 
outflow connection has been inferred from measurements of apparent rotations 
in a number of protostellar jets, although, as noted in Sect. 4.3 below, this 
evidence is still controversial. 

4.3 Centrifugally Driven Winds 

Our focus here is on the role of CDWs as an angular momentum transport 
mechanism. We survey their main characteristics and concentrate on those 

Circular polarization on much larger physical scales was detected along the HH 
135-136 outflow, where it was similarly interpreted as evidence for a helical mag- 
netic field (Chrysostomou et al. 2008). 
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aspects that are relevant to the construction of "combined" disk/wind models, 
the topic of the next subsection. 

Wind Equations 

The qualitative basis for such outflows was considered in Sect. 2.3. In a formal 
treatment, the winds are analyzed using the equations of time-independent, 
ideal MHD.^^ These equations are: 
Mass conservation (continuity equation) 

V • {qv) = , 
Momentum conservation (force equation) 

QvVv = -VP - qV<P + -!- (V X B) X B , (59) 

47r 

where we take the thermal pressure to be given by the perfect gas law 

fimu 

(where T is the temperature, /cb is Boltzmann's constant, fj, is the molecular 
weight and mu is the hydrogen atom mass). 

Induction equation (magnetic field evolution) 

V X {vx B)=0 (61) 

(cf. Eq. (20)), and the solenoidal condition (Eq. (3)). In general, one also needs 
to specify an entropy conservation equation (balance of heating and cooling). 
For simplicity, however, we specialize to isothermal flows (spatially uniform 
isothermal sound speed C = \/P/ g). 

Generic Properties 

Concentrating on axisymmetric flows, it is convenient to decompose the mag- 
netic field into poloidal and azimuthal components, B = Bp + B^, with 
Bp = (Vv4 X 4))/r, where the poloidal flux function A{r, z) can be expressed 
in terms of the poloidal magnetic flux hy A = \I//2'k. The flux function is 
constant along B [B ■ WA = 0) and thus can be used to label the field lines. 
The solution of the induction equation is given by 

X B = Vx , (62) 

The ideal-MHD assumption is justified by the fact that the charged particles' drift 
speeds rapidly become much lower than the bulk speed as the gas accelerates away 
from the disk surface; see Sect. 2.4. 



(58) 



(60) 
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which shows that the electric field E = —{v x B)/c is derivable from an 
electrostatic potential (x)- Since dx/d(j) = on account of the axisymmetry, 
it follows that -Up || Bp. Equivalently, 

QVp = kBp , (63) 

where A: is a flux-surface constant {B ■ Vfc = 0, using V • {gv) = and 
V • -B = 0) . Physically, k is the wind mass load function, 



gvp _ dMvi 



Bp d<I' 



(64) 



whose value is determined by the physical conditions near the top of the disk 
(more precisely, at the sonic, or slow-magnetosonic, critical surface; see discus- 
sion below). By taking the vector dot product of Bp with Eq. (62) and using 
Eq. (63) one finds that the field-line angular velocity Qg = fi ~ (kB^/gr), 
where i7 is the matter angular velocity, is also a flux-surface constant. By 
writing this relation as 

B^ = ^{f2-f2B) (65) 

it is seen that i^s can be identified with the angular velocity of the matter at 
the point where = (the disk midplane, assuming reflection symmetry). 

Using the (j) component of the momentum conservation equation (59) and 
applying Eq. (63) and the field-line constancy of k, one finds that 

l = rv,-'^ (66) 

is a flux-surface constant as well, representing the conserved total specific an- 
gular momentum (matter -|- electromagnetic contributions) along a poloidal 
field line (or streamline). In a CDW the magnetic component of I dominates 
over the matter component near the disk surface, ^ ?'I-Bp^»I ^ ^ whereas 
at large distances this inequality is reversed. The transfer of angular momen- 
tum from the field to the matter is the essence of the centrifugal acceleration 
process. This transfer also embodies the capacity of such a wind to act as an 
eflBcient angular momentum transport mechanism. In fact, as already noted 
in Sect. 2.2, rBz^sB^^s/^i^ represents the magnetic torque per unit area that is 
exerted on each surface of the disk. As was also already discussed in connec- 
tion with magnetic braking, the value of B^^^ is determined by the conditions 
outside the disk, essentially by the inertia of the matter that absorbs the trans- 
ported angular momentum and exerts a back torque on the disk. In the case 
of a CDW, B^s is efi^ectively fixed by the regularity condition at the Alfven 



Eqs. (62)-(65) also hold in a nonideal fluid so long as the field remains frozen 
into a certain particle species (such as the electrons); in this case v and O are 
replaced by the velocity and angular velocity, respectively, of that species (e.g., 
Konigl 1989). 
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critical surface, which is the largest cylindrical radius from which informa- 
tion (carried by Alfvcn waves) can propagate back to the disk. This condition 
yields the value of the "Alfven lever arm" ta, which satisfies 

Z = J7b ri . (67) 

The rate of angular momentum transport by the wind is thus ^ M^Hor^ 
(where we replaced i^s by the midplane value of f2), whereas the rate at which 
angular momentum is advected inward by the accretion disk is Ma.J^o'^o- 
Hence, wind transport can enable accretion at a rate 

Ma^ (rA/ro)'Mw . (68) 

CDW solutions can have rA/^o ~ 3 for reasonable parameters, indicating that 
such outflows could transport the bulk of the disk angular momentum in pro- 
tostellar systems if ~ 0.1 Ma, which is consistent with the observationally 
inferred mass outflow rates. 

By taking the vector dot product of Bp with Eq. (59) , one finds that the 
specific energy 

1 2 tQbB^Bp 
£ = -V — -+W + 9, (69) 

where w = ln(g/gA) (with gA = 47rA:^ being the density at the Alfven 
surface) is the specific enthalpy, is constant along flux surfaces. This is the 
generalized Bernoulli equation, in which the magnetic term arises from the 
poloidal component of the Poynting flux cExB/ An. Using Eqs. (66) and (67) , 
this term can be rewritten as Oq {fior\ — fir'^), and by approximating £ w 
as r oo and assuming that (rA/'^o)^ ^ !> one can estimate the 
value of the asymptotic poloidal speed as 

i^p,co ^ 2i/2r?orA , (70) 

or Wp.co/'i'K ~ 2^/^rA/^'o- This shows that such outflows arc capable of at- 
taining speeds that exceed (by a factor of up to a few) the Keplerian speed 
i;k = r^jfio at their base and thus could in principle account for the measured 
velocities of protostellar jets. 

By combining £, Qb, and I, one can form the field- line constant 

n = £ - QbI=]^v^ -r^^B^ + ^ , (71) 

where we omitted w from the explicit expression. Evaluating at tq gives U = 
— ~%{GM^,QqY/^ , and when this is combined with the form of % at 
large distances (roo) one gets 

rooV4.{r^)Qo - ^-{GM^f'^'nl'^ - ^ [vl{r^) + vl{r^)\ « . (72) 

If one could measure the poloidal and azimuthal speeds at a location Too in 
a protostellar jet and estimate the central mass M*, one would be able to 
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use Eq. (72) (regarded as a cubic in ) to infer the jet launching radius 
ro (GM,//?2)i/3 (Anderson ct al. 2003). One could then also use Eqs. (65) 
and (67) to estimate |-B</,oo /-Spool and rx/ro, respectively, and verify that 
they are significantly greater than 1, as required for self-consistency. This 
has already been attempted in several instances, where it was claimed that 
the results are consistent with a disk-driven outflow that originates on scales 
> 1 AU and carries a significant fraction of the disk angular momentum (e.g., 
Ray ct al. 2007; Coffey ct al. 2007; Chrysostoniou ct al. 2008). However, the 
interpretation of the measurements has been controversial (e.g., Soker 2005; 
Cerqueira et al. 2006) and there have been some conflicting observational 
results (e.g., Cabrit et al. 2006), so the issue is not yet fully settled. 

Critical Points of the Outflow 

The critical points occur in stationary flows at the locations where the fluid 
velocity equals the speed of a backward-propagating disturbance. Since the 
disturbances propagate along the characteristics of the time-dependent equa- 
tions, the critical points can be regarded as relics of the initial conditions in 
a time-dependent flow (Blandford & Payne 1982). If one takes the magnetic 
field configuration as given and considers the poloidal flow along Bp, one can 
derive the location of the critical points and the values of the critical speeds 
(i.e., the values of v-p at the critical points) by regarding the Bernoulli integral 
as a function of the spatial coordinate along the field line and of the den- 
sity ('H = ^{{SjQ), where s is the arc- length of the streamline) and deriving 
the extrema of H by setting dH/dg = (which yields the critical speeds) 
and dH/ds = (which, together with the other relation, yields the critical 
point locations). The critical points obtained from the generalized Bernoulli 
equation (09) correspond to Vp becoming equal to either the slow- or the 
fast-magnetosonic wave speeds (e.g., Sakurai 1985). In the full wind problem, 
the shape of the field lines must be determined as part of the solution by 
solving also the transfield (or Grad-Shafranov) equation, which involves the 
force balance across the flux surfaces. This equation introduces a critical point 
corresponding to Vp = VAp (Okamoto 1975). However, the critical points of 
the combined Bernoulli and transfield equations are in general different from 
those obtained when thcisc two equations are solved separately. The modi- 
fied slow, Alfven, and fast points occur on surfaces that correspond to the 
so-called limiting characteristics (or separatrices; e.g., Tsinganos et al. 1996; 
Bogovalov 1997) . The relevance of the modified critical points has been recog- 
nized already in the original radially self-similar CDW model constructed by 
Blandford & Payne (1982). The modified critical surfaces of the exact, semian- 
alytic wind solutions obtained in this model are defined by the locations where 
the component of the flow velocity that is perpendicular to the directions of 
axisymmetry (i.e., (f>) and self-similarity (i.e., the spherical radius vector .R) 
equals the MHD wave speed in that direction (the 9 direction, using spherical 
coordinates {R, 0, </>}). The significance of the modified fast surface, which in 
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general is located beyond its classical counterpart, is that the poloidal acel- 
eration of the wind continues all the way up to it: initially (roughly until the 
flow reaches the Alfven surface) the acceleration is primarily centrifugal, but 
further out the pressure gradient of the azimuthal field component comes to 
dominate. 

Example: the slow-magnetosonic critical surface 

This critical surface is relevant to the determination of the mass flux in a disk- 
driven wind. There is some uncertainty about this case since the location of the 
first critical point is typically very close to the disk surface (with \z\ being <C 
r), so a priori it is not obvious that ideal MHD is already a good approximation 
there. Under nonideal MHD conditions, all magnetic terms in a disturbance 
are formally wiped out by magnetic diffusivity on its backward propagation 
from "spatial infinity," and one is left with pure sound waves (e.g., Ferreira 
& Pelletier 1995). Here we discuss the situation where the critical point is 
encountered when the charged particles' drift speeds are already small enough 
in comparison with the bulk speed to justify employing ideal MHD. In this case 
the poloidal flow is parallel to the poloidal magnetic field (see Eq. (63)) and 
the relevant wave speed is the slow-magnetosonic (sms) one. We consider the 
standard (rather than the modified) sms point, and our explicit derivation can 
hopefully serve to correct inaccurate statements that have appeared in some 
of the previous discussions of this topic in the literature. 

In the magnetically dominated region above the disk surface the shape of 
the field lines changes on the scale of the spherical radius R. Anticipating that 
the height ^sms of the sms point is <C R, we approximate the shape of the 
field line just above the point {vg, Zs} on the disk surface by a straight line 
(cf. Sect. 2.3): 

r = + s sin 6^ , z = + s cos 6g , (73) 

where the angle 0s gives the field-line inclination at the disk surface (sin 9s = 
Br,s/Bpy, tan^s = Br,s/-Bz)- The Bernoulli integral (in the form of Eq. (71)) 
then becomes, after substituting = -GM*/(r^ + z^YI"^, 



+ C2ln(^), (74) 



{(ts + ssin6's)2 + (zg + s cos 6*8)2] 1/2 
where f2{s, g) is given by combining 



1 - 



^ = (75) 
a 

(obtained from Eqs. (65) (67)) and Eq. (73). In the hydrostatic approximation 
to the disk structure (see Sect. 4.4) 1?b = 17K(rs). We expect this equality to 
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hold approximately also for the exact solution, in which Vz 0, so that, in 
particular, » i7 in the region of interest. 

Setting cyU/dg = yields the speed of an sms wave propagating along Bp: 

Vs^, = 4F C' , (76) 

where = (-B^ ,.+5^ Ji/2 and we approximated {B^^sms/ B^msf by [B^^s/Bsf 
and assumed q/ qa ^ 1 and {Q/Qb — ^Y ~ (^A/7')^(£'A/^?)^ ^ 1 in the region 
between the top of the disk and the sms surface. Setting also dH/dr = and 
approximating rgms ~ ''s then gives the height of this point: 

'2'sms 3 tan (77) 



3 tan^ 9.-1 



Equation (77) yields a meaningful result only if ianOs > l/\/3; i-c., if the field 
line is inclined at an angle > 30° to the z axis. This is the CDW launching 
condition in a Keplerian disk (Eq. (6)) that was derived in Sect. 2.3 using the 
mechanical analogy to a "bead on a rigid wire." The relationship between these 
two derivations becomes clear when one notes that, in the limit {f2 / f^B — ^)^ «C 
1, the 2nd term on the right-hand side of Eq. (74) becomes equal to the 
centrifugal potential — i7^r^/2, so that the 2nd and 3rd terms (which together 
dominate the right-hand side of this equation) are just the effective potential 
^eff used in Eq. (5).^^ The correspondence of the sms point to the maximum 
of <?eff was already noted in Sect. 2.3. 

The density at the sms point can be related to the density at the disk 
surface by evaluating the energy integral (74) at both Zg and ^;sms- This yields 

-^=''^[-2-2[-^) 3tan^g.-l |' ^'^^ 

where we assumed (wp.s/i'sms)^ 1. The mass flux injected into the wind 
from the two sides of the disk is then 

'^QsmsVsmsCOsOs = 2QsmsC-^ (79) 



n 1 — trains^ aula —traiiia^ 7-> 

27rr dr Bg 

(e.g., Lovelace et al. 1995). 

When 9s approaches (and decreases below) 30° the field-line curvature 
needs to be taken into account in the analysis. The height of the sonic point 
rapidly increases to ~ R, with the potential difference growing to ~ GM^,/R: 
the launching problem becomes essentially that of a thermally driven spherical 
wind (e.g., Levinson 2006). 



Furthermore, the condition dH/dr = is equivalent to Bp- V<feff = in this limit 
(Campbell 2002), so the latter relation can also be used to obtain the location of 
the sms point in this case (e.g., Ogilvie 1997). 
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It has been argued (Spruit 1996) that the field-Hne inclination increases 
systematically with radius along the disk surface and that only in a narrow 
radial range the conditions are favorable for driving a wind that both satisfies 
the CDW launching condition and is not "overloaded" (and hence conceiv- 
ably highly unstable; cf. Cao & Spruit 1994). This is an intriguing suggestion, 
given the fact that so far there is no observational evidence for an extended 
wind-driving region in protostellar disks (see Sect. 4.1). However, in principle 
it may be possible to launch oiitflows over a large radial range. In particular, 
as illustrated by the similarity solution presented in Sect. 3, if the magnetic 
flux is advected inward by the accretion flow, the field-line inclination could be 
favorable for launching and need not change strongly along the disk. Further- 
more, evidently the mass loading of stable disk/wind configurations actually 
decreases as 6s is increased (see Sect. 4.5). 



Exact Wind Solutions 



As was already noted in Sect. 2.5, the radial self-similarity approach has been 
used to construct exact global solutions of CDWs. The basic character of 
this model is revealed by the prototypal Blandford & Payne (1982) solution. 
The underlying assumption that all quantities scale as a power law of the 
spherical radius implies that all the critical surfaces are conical. Furthermore, 
all the relevant speeds (including the fluid, sound, and Alfven) must scale like 
the characteristic speed of the problem {vk), i-e., as i?^^/^. The scaling of 
the magnetic field amplitude can be inferred from the vertically integrated 
thin-disk equations presented in Sect. 3.2, in which we now set d/dt = 0. In 
particular, if the mass outflow from the disk has only a negligible effect on 
the accretion rate then Ma = 27rr|i;r|^ — const by Eq. (31) and the angular 
momentum conservation equation (33) in a Keplerian disk can be written as 

^M^VK = r^\B,B^^s\ . (80) 

Equation (80) implies the similarity scaling B oc r~®/^, from which we infer 
by dimensional arguments (wa oc B j ^Jq) that q cx r~^/^ in this case. It then 
follows that Afw oc Inr (cf. Eq. (79)), which is consistent with the mider lying 
assumption that only a small fraction of the inflowing mass leaves the disk 
over each decade in radius. More generally, one can define a mass ejection 
index ^ > by 

dlnr 

(e.g., Ferreira & Pelletier 1995) and deduce 

ding ^e_5 dln^^ 3 

dlnr 2 4 ' dlnr ^2 ^ ^ 
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(e.g., Contopoulos & Lovelace 1994), with ^ — >■ corresponding to the 
Blandford & Payne (1982) solution. The poloidal electric current scales as 
/ oc rB^ oc r^^i~^')/'^. For ^ > 1/2 the flow is in the current-carrying regime, 
with the poloidal current density being antiparallel to the magnetic field. In 
this case the current tends to zero as the symmetry axis is approached, so such 
solutions should provide a good representation of the conditions near the axis 
of a highly coUimated flow. Conversely, solutions with ^ < 1/2 correspond to 
the return-current regime (in which the poloidal current density is parallel to 
the field) and are most suitable at larger cylindrical distances. The Blandford 
& Payne (1982) solution is sometimes critiqued for having a singular behavior 
on the axis."'^* One should bear in mind, however, that even though the de- 
tailed global current distribution (which includes both a current-carrying and 
return-current regimes) cannot be represented by the simplified self-similar 
solution, this flaw is not fundamental. Besides, given that the disk has a flnite 
inner radius, the issue of the behavior of the flow at r = is, to a certain de- 
gree, merely academic. Whereas the value of the ejection index is arbitrary in 
pure wind models, it becomes an eigenvalue of the self-similar solution (fixed 
by the regularity condition at the sonic, or sms, critical surface; see Sect. 4.3) 
when one considers a "combined" disk/wind model (e.g., Li 1996a). 

While the magnetic field lines must bend away from the symmetry axis 
within the disk in order to satisfy the wind launching condition (6) at the 
surface, as soon as the magnetically dominated region at the base of the wind 
is reached thej^ start to bend back toward the axis on account of the magnetic 
tension force (see Sect. 2.1). Further coUimation is achieved in current-carrying 
jets by the hoop stress of the toroidal magnetic field (the term —J^B^/c = 
— {l/^Tir'^)d{rB^Y/dr in the radial force equation), the analog of a z-pinch in 
laboratory plasma experiments. The asymptotic behavior in general depends 
on the current distribution: if / — )• as r — )• oo then the field lines are space- 
filling paraboloids, whereas if this limit for the current is finite then the fiow is 
coUimated to cylinders (e.g., Heyvaerts & Norman 1989). In practice, however, 
self-similar wind solutions typically do not reach the asymptotic regime but 
instead sc^lf-focus (with streamlines intersecting on the axis) and terminate at 
a finite height (e.g., Vlahakis et al. 2000). 

The cold self-similar wind solutions are specified by any two of the follow- 
ing three parameters: 



Interestingly, one can also generalize the Blandford & Payne (1982) solution to a 
class of semianalytic, but non-self-similar solutions, defined by the constancy of a 
certain function of the magnetic flux that controls the shape of the Alfven surface 
(PcUoticr & Pudritz 1992). The Blandford & Payne solution then separates wind 
configurations in which the initial field-line inclination increases progressively with 
radius from outflows that emerge from a bounded region and in which the initial 
field inclination decreases with r and the field lines converge into a cylindrical 
sheath. 

It is in fact a double singularity, since on the axis itself there is an oppositely 
directed line current that exactly compensates for the distributed return current. 
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A 



normalized specific angular momentum 



K 




normalized mass/magnetic flux ratio 



(83) 



6r,s = ^5— poloidal field inclination at the disk surface , 

where br^s must satisfy the constraint (6). A viable solution is further charac- 
terized by A > 1 (typically » 1) and k < 1 (typically <C 1). A two-parameter 
choice yields a solution if the corresponding flow crosses the modifled Alfven 
critical surface (see Sect. 4.3); if the two parameters are k and A then, for any 
given value of k, this can happen only if A exceeds a minimum value that ap- 
proximately satisfies KAmin(A,„in — 3)"^/^ = 1.-^^ In Sect. 4.4 we discuss joining 
a global wind solution of this type to a radially localized wind-driving disk 
solution. In this case the sonic (or sms) critical-point constraint imposed on 
the latter solution yields the value of k, and, in turn, the Alfven critical-point 
constraint on the wind solution fixes one of the disk model parameters (for 
example, by providing the value of B^^s/Bz,s = — '«(A — 1) from given values 
of K and 6r,s)-^'^ 

Observational Implications 

The exact wind solutions and estimates of the physical conditions around 
young stellar objects make it possible to identify a number of physical char- 
acteristics of disk-driven protostcUar winds that could have potentially sig- 
nificant observational implications. While this topic is not directly within the 
scope of this book, we nevertheless briefly describe some of these properties 
here inasmuch as it may not always be feasible in practice to distinguish be- 
tween the observational signatures of a disk and a wind. Furthermore, invoking 
a disk-driven wind may sometimes provide a straightforward explanation to a 
disk observation that may otherwise seem puzzling. Besides the references ex- 
plicitly listed below, the reader may also consult the review articles by Konigl 
& Ruden (1993) and Konigl & Pudritz (2000) on this subject. 

Centrifugal driving is found to be an efficient acceleration mechanism, 
leading to a rapid increase in the poloidal speed, and a correspondingly strong 
decrease in the wind density, above the disk surface (see Fig. 4).^^ This leads 

This expression differs slightly from the one that appeared in Blandford & Payne 
(1982). 

Although self-similar CDW solutions in which the wind passes also through the 
fast-magnetosonic critical surface have been obtained, the implications of this 
additional constraint to a global disk/wind model remain unclear. Ferreira & 
Casse (2004), for instance, suggested that in order to cross this surface the outflow 
must experience significant heating after it leaves the disk and that the added 
constraint is related to this requirement. 

The model wind presented in Fig. 4 exhibits a comparatively weak collimation. 
The degree of collimation depends on the model parameters and is particularly 
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to a strong stratification, which, in fact, is most pronounced in the z direction 
even in the i?-self-siniilar wind models (Safier 1993b). The strong momentum 
flux in the wind also implies that an outflow originating from beyond the 
dust sublimation radius could readily uplift dust from the disk (Safier 1993a). 
This may load to viewing anglc-dcpcndcnt obscuration and shielding of the 
central continuum , whose effect would be distinct from that of molecular cloud 
gas undergoing gravitational collapse. The stratified dust distribution could 
intercept a large portion of the central continuum radiation and reprocess 
it to the infrared, which could be relevant to the interpretation of the IR 
spectra of low- and intermediate-mass star/disk systems (e.g., Konigl 1996; 
Tambovtscva & Grinin 2008). Furthermore, scattering by the uplifted dust 
could contribute to the observed polarization pattern in these objects. 

Although, as we already noted, ideal MHD should be an excellent approxi- 
mation for modeling the dynamics of the flow (except perhaps right above the 
disk surface, where the nonideal formalism employed within the disk might 
still be relevant for determining the first critical point of the wind), the esti- 
mated degree of ionization at the base of a protostellar disk wind is typically 
low enough that energy dissipation induced by ion-neutral drag (ambipolar 
diffusion) could play an important role in the wind thermodynamics. This is 
because the volumetric heating rate scales as |(V x B) x Bp/7iPip oc 
(see Sect. 2.4), where the distributions of B and g are regarded as being 
fixed by the flow dynamics, whereas adiabatic cooling, the most important 
temperature-lowering mechanism, remains relatively inefficient in a coUimated 
disk outflow. Under these circumstances, the temperature might rise rapidly 
(Safler 1993a), although its precise terminal level depends on the relevant 
value of the coUisional drag coefficient 71 (Shang ct al. 2002). This heating 
mechanism could contribute to the observed forbidden line emission (Safier 
1993b; Cabrit et al. 1999; Garcia et al. 2001a,b) and thermal radio emission 
(Martin 1996) from CTTSs, although in both of these applications an addi- 
tional source of heating might be needed, possibly associated with dissipation 
of weak shocks and turbulence (e.g., O'Brien et al. 2003; Shang et al. 2004). 

On large scales, centrifugally driven outflows assume the structure of a 
coUimated jet (most noticeable in the density contours) and a surrounding 
wide-angle wind (Shu et al. 1995; Li 1996b). This bears directly on the general 
morphology of protostellar sources, and in particular on the shapes of the 
(wind/jet)-driven outflow lobes that form in the surrounding molecular gas. 

sensitive to the mass loading of the outflow (e.g., Pudritz et al. 2007). Protostellar 
jets are inferred to correspond to more highly coUimated disk winds, in which the 
enhanced mass loading might be attributed to the presence of a disk corona (e.g., 
Dougados et al. 2004). 

In an application to active galactic nuclei, this effect was invoked by Konigl & 
Kartje (1994) to account for the obscuring/absorbing "molecular torus" identified 
in the centers of Seyfert galaxies. 

For an application of this model to the polarization properties of Seyfert galaxies, 
see Kartje (1995). 
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r(AU) 

Fig. 4. Structure of a Blaridford-Payne-type self-similar CDW from a protostellar 
disk (model C in Safier 1993a). Top: Meridional projections of the streamlines (la- 
beled by the value, in AU, of the radius where they intersect the disk). Bottom: 
Contours of Up,oo/'^Ki {solid lines), where Vp^ao is the asymptotic poloidal speed and 
vki is the Keplerian speed at 1 AU, and of log (g/ qi) {dashed lines), where qi is the 
density at the base of the wind at 1 AU. The heavy line on the left indicates the 
streamline that originates at 0.1 AU. 

4.4 Equilibrium Disk/Wind Models 

In contrast with the wind zone, where the dynamics is well described by 
an effectively infinite conductivity, the accretion disk in protostellar systems 
must be modeled using nonideal MHD. The structure of the accretion flow 
thus depends critically on the properties of the conductivity tensor at each 
point, which are in turn determined by the spatial distribution of the degree of 
ionization (see Eqs. (12)-(14)). We therefore start this subsection by reviewing 
the physical processes that affect the disk ionization (see also Chap. II [The 
Chemical Evolution of Protoplanetary Disks] ) , which would allow us to choose 
the relevant physical parameters for our models. In reality the accretion flow 
itself has an influence on the ionization structure — through its effect on 
the disk column density (see Sect. 2.4) or on the distribution of dust grains 
(see Sect. 6), for example — which a truly self-consistent model must take 
into account. So far only simpler models have been constructed, in which 
this influence is not fully accounted for and a variety of approximations are 
employed. We discuss these models next, describing the basic properties of the 
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derived equilibrium solutions and distinguishing between "strongly coupled" 
and "weakly coupled" disk configurations. 

Disk Ionization Structure 

The dependence of the degree of field-matter coupling on the abundances of 
the ionized species in the disk (ions, electrons and charged dust grains) is a 
direct outcome of the fact that this coupling is effected by the collisions of these 
particles with the much more abundant neutrals. The degree of ionization is 
given by the ratio of the sum of the positively (or, equivalently, negatively) 
charged particle densities to the neutral particle density. It is calculated by 
balancing the ionization and recombination processes operating in the disk. 

Interstellar (and possibly also protostellar) cosmic rays as well as X-ray and 
far-UV radiation produced by the magnetically active protostar (e.g.. Glass- 
gold et al. 2000, 2005) are the main potential ionizing agents in protostellar 
environments. It is, however, unclear how effective cosmic-ray ionization really 
is, since the low-energy particles most relevant for this purpose may be de- 
flected by magnetized disk or stellar outflows or by magnetic mirroring near 
the disk surface, and they might also be scattered by magnetic turbulence 
within the disk. If present, they may be more important than X-rays in the 
inner disk (r < 1 AU) , where the surface density is larger than the X-ray at- 
tenuation length. But if cosmic rays are excluded from this region and ionized 
particles are not transported there by othcir means (e.g., Tiu-ner et al. 2007), 
the gas near the midplane might be ionized only at the low rate of decay of 
radioactive elements such as ^°K. On the other hand, free electrons are rapidly 
lost through recombination processes that occur both in the gas phase (where 
electrons recombine with molecular and metal ions via dissociative and radia- 
tive mechanisms, respectively) and on grain surfaces, and through sticking to 
dust particles (e.g., Oppenheimer & Dalgarno 1974; Nishi et al. 1991). The 
ionization equilibrium is, however, very sensitive to the abundance of metal 
atoms, as they rapidly remove charges from molecular ions but then recom- 
bine much more slowly with electrons (e.g., Sano et al. 2000; Promang et al. 
2002). 

Dust grains can also significantly affect the degree of field-matter coupling 
if they are mixed with the gas. They do so in two ways. First, they reduce the 
ionization fraction by absorbing charges from the gas and by providing addi- 
tional recombination pathways for ions and electrons. Second, in high-density 
regions the grains themselves can become an important charged species (e.g., 
Nishi et al. 1991), leading to a reduction in the magnetic coupling because 
grains have much smaller Hall parameters than the much less massive elec- 
trons and ions. Wc note in this connection that, if dust grains settle; to the 
midplane, the ionization fraction may become sufficiently large to provide ad- 
equate magnetic coupling even in the inner (< lAU) disk regions (Wardle 
2007). This conclusion could, however, be mitigated in the prcisence of turbu- 
lence, which might leave a residual population of small dust grains (carrying 
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a significant fraction of the total grain charge) suspended in the disk (e.g., 
Nomura & Nakagawa 2006; Natta et al. 2007). 

Exact Disk Solutions 

To gain physical insight into the distinguishing properties of wind-driving 
disks, we consider a simplified model originally formulated by Wardle & Konigl 
(1993). In this model, the entire angular momentum of the accreted matter is 
assumed to be transported by a CDW. The disk is taken to be geometrically 
thin, vertically isothermal, stellar gravity-dominated, and threaded by an open 
magnetic field (possessing reflection symmetry about the midplane). Under 
the thin-disk approximation, the vertical magnetic field component is taken 
to be uniform with height in the disk solution. The disk gas is assumed to be in 
the ambipolar diffusivity regime, with ions and electrons being the dominant 
charge carriers. In the original model, the ion density was taken to be constant 
with height (which could be a reasonable approximation in the outer regions 
of certain real systems). The main simplifications involve considering only 
a radially localized (Ar <C r) disk region and (in accord with the thin-disk 
approximation) retaining the z derivatives but neglecting all r derivatives 
except those of v^f, (which scales as r~^/^) and Br (which appears in V • B = 
0). The neglect of the radial derivatives in the mass conservation equation 
(Eq. (58)) implies that the vertical mass flux gvz is taken to be uniform 
with height in the disk solution. This solution is extended through the first 
critical point (which is the thermal sonic point if the flow is still diffusive and 
the sms point if it is already in the ideal-MHD regime; see Sect. 4.3) and then 
matched onto a global Blandford & Payne (1982) ideal-MHD wind solution. 
These simplifications do not compromise the physical essence of the results. In 
particular, the qualitative characteristics remain unchanged when the radially 
localized disk solution is generalized to a global configuration in which both 
the disk and the wind are described by a single self-similar model that includes 
both z and r derivatives (e.g., Li 1996a). Furthermore, solutions with similar 
properties are obtained when the disk is in the Hall or Ohm diffusivity regimes 
(e.g., Salmeron et al. 2010) and when the full conductivity tensor is used in 
conjunction with a realistic ionization profile (see Fig. 5). 

This approximation has been critiqued (e.g., Ferreira 1997) for not allowing Vz 
to assume negative values within the disk, as it must do in cases (expected to 
be typical) in which the disk thickness decreases as the protostar is approached. 
This issue can be fully addressed only in the context of a global disk/wind model. 
However, it can be expected that any error introduced by this approximation 
would be minimized if the upward mass flux remained small enough for Vz to 
have only a weak effect on the behavior of the other variables within the disk. 
Under the assumption that \vr\ is of the same order of magnitude as — vk\ 
one can readily show that the condition for this to hold is that v^/C remain <C 1 
everywhere within the disk. This can be checked a posteriori for each derived 
solution. 
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Fig. 5. Vertical structure of a radially localized wind-driving disk solution at 1 
AU around a Sun-like protostar. The disk is assumed to have a column density of 

S = 600gcm~^ and its model parameters arc ao = 0.75, C/vk = 0.1, e — 0.1, and 
£b = 0. The parameters of the matched radially self-similar CDW are k = 2.6 x 10~^, 
A = 4.4 X 10^, and br,s = 1.6. (Left) Top: Ambipolar, Hall, and Ohm diffusivities. 
Bottom: Ionization rates by cosmic rays (cr), X-rays (xr), and radioactivity (rad), 
and electron fraction rie/riH. (Right) Top: Radial and azimuthal magnetic field 
components, mass density, and Elsasser number. Bottom: Velocity components. 
The mass accretion rate for this model is 7 x 10~^ Mq yr~^, which is consistent with 
the inferred values for the early (Class 0/Class I) protostellar accretion phase. 



In general, an equilibrium disk solution is specified by the following pa- 
rameters: 

1. ao = VAo/C, the midplane ratio of the Alfven speed (based on the uniform 
vertical field component) to the sound speed. This parameter measures the 
magnetic field strength.^''' 

2. C/vk, which in a thin isothermal disk is equal to hT/r, the ratio of the 
tidal (i.e., reflecting the vertical gravitational compression) density scale 
height to the disk radius. While this parameter, which measures the geo- 
metric thinness of the disk, does not appear explicitly in the normalized 
structure equations, it nevertheless serves to constrain physically viable 
solutions (sec Eq. (84) below). 

3. The midplane ratios of the conductivity tensor components: [iJp/c7_l]o (or 
[c7h/(t_l]o) and [(t_l/(To]o- They characterize the conductivity regime of 
the gas (see Sect. 2.4). When there are only two charged species (ions 

Note that ao is related to the midplane plasma beta parameter /3o through oo = 
(2/7/3o)^'^^, where 7 denotes the adiabatic index of the fluid. 
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and electrons) , one can equivalently specify the midplane Hall parameters 

/3io and /3co- In the inner disk regions the conductivity tensor components 
typically vary with height, reflecting the ionization structure of the disk 
(see Fig. 5). 

4. The midplane Elsasscr number Aq (sec Eqs. (23) and (25)), which mea- 
sures the degree of coupling between the neutrals and the magnetic field. 

5. e = —Vrfi/C, the normalized inward radial speed at the midplane. Al- 
though the value of e could in principle be negative (as, in fact, it is in 
certain viscous disk models; e.g., Takeuchi & Lin 2002), it is expected 
to remain > when a large-scale magnetic field dominates the angular 
momentum transport. 

6. cb = —VBrfl/C, the normalized radial drift speed of the poloidal magnetic 
field lines (see Eq. (22)). This parameter vanishes in a strictly steady-state 
solution but is nonzero if, as expected, the magnetic field lines drift radially 
on the "long" accretion time scale r/\vr\- It is incorporated into the model 
through the z component of Eq. (20) (i.e., one sets dB^/dt ^ but keeps 
dBr/dt = and dB^/dt = 0)?^ 

We already remarked in Sect. 4.3 on how the wind solution, through the 
Alfven-surface regularity condition, can be used to constrain the disk solution 
when the two are matched. In particular, if the disk solution is used to fix 
the wind parameters k (through the sonic/sms critical-point constraint) and 
6r,sj the wind solution yields B^^i^/B^, which in turn can be used to obtain 
the value of the disk parameter e. Previous "combined" disk/wind models 
treated ee as a free parameter and typically set it equal to zcro.^^ Under this 
approach, 6^,8 was fully determined by the conditions inside the disk. However, 
as we noted in Sect. 3.2 (see also Ogilvie & Livio 2001), br^s is in fact also 
determined by the conditions outside the disk and can be directly related to 
the distribution of B;^ along the disk if the force-free field above the surface 
can be adequately approximated as being also current free (see Eq. (40)). In 



It is worth noting, though, that the poloidal components of Eq. (20) can be written 
as d'P/dt = -2nrBoVBrA), with = {l/2-Kr)d<P/dr and Br = -{l/2-Kr)d^ /dz. 
Thus, if Bz changes because of the slow radial diffusion of the flux surfaces, so 
will also Br- However, it can be argued (Konigl ct al. 2010) that Br (and similarly 
also B^) can in principle change at any given radial location on the much shorter 
dynamical time r/v^ through field-line shearing by the local velocity field. One 
can therefore assume that Br and B^ attain their equilibrium configurations on 
the "short" time scale and neglect, in comparison, the much slower variations that 
are expressed by the explicit time derivate terms. The B^ component is distinct 
in that it can only change on the "long" radial drift time. 

Wardle & Konigl (1993) noted, however, that for physical consistency one has to 
require es < e. They also demonstrated that the solution variables (except for 
B^^a and Vr,o) are insensitive to the value of this parameter, as expected from 
the fact that the only modification to the equations introduced by varying tB 
involves changing the radial velocity of the reference frame in which the poloidal 
field lines are stationary. 
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a more general treatment of the disk/wind problem, this constraint can be 
used to fix the value of (soe Teitlcr 2010). 

Before turning to the specific solution displayed in Fig. 5, it is instructive 
to list the dominant magnetic terms (under the thin-disk approximation) in 
the neutrals' force equation (Eq. (59)) and review their main effects inside the 
disk. We have (sec Eqs. (28)-(30)): 

• radia/ component: 

representing the magnetic tension force that acts in opposition to central 

gravity; 

• azimuthal component: ^-g^ 

representing the magnetic torque that transfers angular momentum from 

the matter to the field; 

• wemcai component: — ^ — stt 

representing the magnetic squeezing of the disk (which acts in the same 
direction as the gravitational tidal force and in opposition to the thermal 

pressure-gradient force) . 

The solution shown in Fig. 5 describes the vertical structure of a disk that 
changes from being Hall-dominated near the midplane to being ambipolar- 
dominated near the surface. While the qualitative properties of the solution 
are not sensitive to the details of the diffusivity profile, it is heuristically 
useful to consider a distinct diffusivity regime. We therefore specialize in the 
following discussion to the "pure ambipolar" (with gi{z) = const) case treated 
in Wardle & Konigl (1993). As illustrated in Fig. 6 (in which it is assumed 
that eB — 0), the disk can be vertically divided into three distinct zones:^® a 
quasi-hydrostatic region near the midplane, where the bulk of the matter is 
concentrated and most of the field-line bending takes place, a transition zone 
where the inflow gradually diminishes with height, and an outflow region that 
corresponds to the base of the wind. The first two regions are characterized 
by a radial infiow and sub-Keplerian rotation, whereas the gas at the base of 
the wind flows out with > Vk- 

• The quasi-hydrostatic region is matter dominated, with the ionized parti- 
cles and magnetic fleld being carried around by the neutral material. The 
ions are braked by a magnetic torque, which is transmitted to the neu- 
tral gas through the frictional drag; therefore Vi^ < in this region. The 
neutrals thus lose angular momentum to the field, and their bac;k reaction 
leads to a buildup of the azimuthal field component away from the mid- 
plane. The loss of angular momentum enables the neutrals to drift toward 
the center, and in doing so they exert a radial drag on the field lines. This 
drag must be balanced by magnetic tension, so the field lines bend away 
from the rotation axis. This bending builds up the ratio Br/B^, which 

Our discussion pertains to disks in which yl > 1 throughout their entire vertical 
extent, but is also applies to the A > 1 wind-driving surface layers of "weakly 
coupled" disks. 
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Fig. 6. Schematic diagram of the vertical structure of an ambipolar diffusion- 
dominated disk, showing a representative field line and the poloidal velocities of 
the neutral and the ionized gas components. Note that the poloidal velocity of the 

ions vanishes at the midplanc (z = 0), consistent with the assumption that es ~ 0, 
and that it is small for both components at the top of the disk (z = Zs)- The relation 
between the azimuthal speeds is also indicated. 



needs to exceed l/\/3 at the disk surface to launch a centrifugally driven 
wind. The magnetic tension force, transmitted through ion-neutral colli- 
sions, contributes to the radial support of the neutral gas and causes it to 
rotate at sub-Kcplcrian speeds. 

• The growth of the radial and azimuthal field components on moving away 
from the midplane results in a magnetic pressure gradient that tends to 
compress the disk. The magnetic energy density comes to dominate the 
thermal and gravitational energy densities as the gas density decreases, 
marking the beginning of the transition zone (at z h, where h is the 
density scale height). The field above this point is nearly force free and 
locally straight (see Sects. 2.1 and 2.3). 

• The field angular velocity /?b is a flux-surface constant (see Eq. (65) and 
footnote 15; note, however, that this strictly holds only when cb is iden- 
tically zero). The ion angular velocity Vi^/r differs somewhat from f2B 
but still changes only slightly along the field. Since the field lines bend 
away from the symmetry axis, the cylindrical radius r, and hence Vi^, 
increase along any given field line, whereas i;^ decreases because of the 
near-Keplerian rotation law. Eventually a point is reached where (vi^ — v^) 
changes sign. At this point the magnetic stresses on the neutral gas are 
small and its angular velocity is almost exactly Keplerian.-^^ Above this 



This result was utilized in the derivation of the slow-magnetosonic critical surface 
properties in Sect. 4.3. 
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point the field lines overtake the neutrals and transfer angular momentum 

back to the matter, and the ions start to push the neutrals out in both the 
radial and the vertical directions. This region can be regarded as the base 
of the wind, and one can accordingly identify the disk surface with the 
location where f0 becomes equal to wk- The mass outflow rate is fixed by 
the density at the sonic/sms point (marked by Zsms on the diagram). 

The Hydrostatic Approximation 

The structure equations of the radially localized disk model can be simplified 
by setting ~ 0. Although this approximation is most appropriate for the 
quasi-hydrostatic region, one can nevertheless extend it to the disk surface 
to obtain useful algebraic constraints on viable disk solutions. For the pure 
ambipolar regime (and again assuming gi = const), they are given by (Wardle 
& Konigl 1993; Konigl 1997) 

(2ro)-i/2 < ao < V3 < eTo < vk/2C , (84) 

where the parameter Tq represents the midplane Elsassser number in the 
ambipolar limit (see Eq. (25)). The four inequalities in Eq. (84) have the 

following physical meaning (going from left to right): 

1. The disk remains sub-Kcplcrian everywhere below its surface. 

2. The wind launching condition (6^,3 > 1/v^) is satisfied. 

3. The top of the disk (zg) exceeds a density scale height {h « (ao/eTg) h^), 
ensuring that the bulk of the disk material is nearly hydrostatic and that 
Mw remains < Ma.^° 

4. The midplane Joule heating rate (see Sect. 2.4) is less than the rate of 
gravitational potential energy release there {Qo\vrfi\v^/'^f per unit vol- 
ume), as the latter is the ultimate source of energy production in the 
disk. 

The inequalities (1) and (2) together place a lower bound on the neutral- 
field coupling parameter Tq, for which a more detailed analysis yields the 
value of ~ 1 (see Konigl et al. 2010). The inequalities (2) and (3), in turn, 
imply that h/hq^ < 1, i.e., that magnetic squeezing dominates tidal gravity. 
Analogous constraints on wind-driving disks can be obtained also in the Hall 
and Ohm diffusivity regimes. For example, in the Hall case the parameter 
space divides into four sub-regimes with distinct sets of constraints (Konigl et 
al. 2010; Salmeron et al. 2010). In all of these cases it is inferred that Tq > 1 
and /i//iT < Ij which indicates that these are generic properties of wind- 
driving disks. In practice one finds that successful full solutions typically have 
ao < 1 (the midplane magnetic pressure is smaller than the thermal pressure, 
but not greatly so) and e < 1 (the midplane inflow speed is not much smaller 
than the speed of sound). 

These requirements also place upper limits on the ratio of the density at the 
sonic/sms point to the midplane density (see Eq. (78) for the sms case). 
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The relations listed in Eq. (84) have a few other interesting implications. 
For example, one can use the inequality (3) to argue that wind-driving disks 
are stable to the fastest growing linear mode of the MRI, because it implies 
that the vertical wavelength of this mode, ~ vao/^k (see Chap. V), is larger 
than the magnetically reduced disk scale height h = (ao/eTo)C/i7K- But one 
can also argue that the inequality (1) leads to a useful criterion for the onset 
of the MRI in the disk. Such a criterion is of interest since, as we already 
noted when the Elsasscr number was introduced in Sect. 2.4, the minimum- 
coupling condition for the development of MRI turbulence in a diffusive disk 
(represented by a lower bound on A) is evidently similar to that inferred for 
driving a disk wind. The question then arises as to whether, in a disk that is 
threaded by a large-scale, ordered field, both vertical (wind-related) and ra- 
dial (MRI-related) angular momentum transport can occur at the same radial 
location; this can alternatively be phrased as a question about the maximum 
magnetic field strength for the operation of the MRI. When the inequality 
(1) is violated, the surface layers become super-Keplerian, implying outward 
streaming motion that is unphysical in the context of a pure wind-driving 
disk. However, as elaborated in Salmeron et al. (2007), such motion could be 
associated with the two-channel MRI mode that underlies MRI-induced tur- 
bulence. If one regards the parameter combination 2Ta^ that figures in this 
inequality as a function of height in the disk rather than being evaluated at 
the midplane, one can infer from the fact that it scales as B^{gi/g) that it 
generally increases with z. It is thus conceivable that the inequality (1) (gen- 
eralized in this manner) is violated near the midplane but is satisfied closer to 
the disk surface. Salmeron et al. (2007) developed disk models in which both 
radial turbulent transport and vertical transport associated with the m,ean 
field take place in the region where the inequality (1) is violated but only 
vertical transport (with ultimate deposition of the removed angular momen- 
tum in a disk wind) occurs at greater heights. They concluded, however, that 
significant radial overlap between these two transport mechanisms is unlikely 
to occur in real disks. 

Weakly Coupled Disks 

In the disk models considered so far, the minimum-coupling condition on the 
neutrals, yl > 1, was satisfied throughout the vertical extent of the disk. 
We refer to such disks as being strongly coupled. As we have, however, just 
noted, the parameter values in a real disk could vary with height, reflecting 
the vertical stratification of the column density (which shields the ionizing 
radiation or cosmic rays) and the density. In particular, if the disk is in the 
ambipolar regime near the surface, in the Hall regime further down, and in 
the Ohm regime near the midplane, then A scales as Qi, Qi/g, and gi/g"^, 
respectively, on going from z ~ Zs to z = (see Eq. (25)). The Elsasser 
number will thus increase with height on moving up from z = as the gas 
becomes progressively more ionized and the density decreases. It will generally 
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peak on reaching the ambipolar regime and will subsequently drop as Qi (which 
typically scales as £i to a power between and 0.5) decreases. The gas will be 
weakly coupled in regions where A remains <C 1- Weakly coupled disks have 
distinct properties from those of strongly coupled ones, as first pointed out 

by Li (1996a). 

As an example, consider a protostellar disk that is ionized solely by cosmic 
rays (with the ionization rate decreasing exponentially with depth into the 
disk with a characteristic attenuation column of 96gcm^^) and in which the 
charge carriers are small, singly charged grains of equal mass (a reasonable 
approximation for the high-density inner regions of a real disk; e.g., Neufeld 
& Hollenbach 1994). Near the disk surface the gas is in the ambipolar regime, 
but if the disk half-column is :» 96gcm~^, the degree of ionization near the 
midplane is very low and the gas there is in the Ohm regime. Figure 7, taken 
from Wardle (1997), depicts two illustrative solutions obtained for different 
radii in this model. The solution on the left corresponds to a large enough 
radius (and, correspondingly, a sufficiently low disk column) for the disk to 
be in the ambipolar regime throughout its vertical extent. By contrast, the 
solution on the right depicts a smaller radius where the column density is large 
and the disk is weakly ionized (and in the Ohm regime) near the midplane. 



0.5 



.\2p/po 










-B(|)/Bz ' 










: Z 










V-Vk$ 

C : 



1 1 1 1 


1 1 1 1 1 1 1 . 




' -B^Bz ""^ 


\ 4p/po 


/ 




/ 




/ B,/Bz " 




/ 














/ Z J 






V-V|($ 




C 

1 





1 2 
z/hj 



3 



2 4 
z/hj 



Fig. 7. Strongly coupled (left) vs. weakly coupled (right j disk solutions. 



The main differences between these two solutions can be summarized as 
follows: 



Owing to the assumed equal mass of the positive and negative charge carriers, 
the Hall term in Ohm's Law is identically zero (see Eq. (13)). 
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Strongly coupled disks 

• ^^Ao ^ C (midplane magnetic pressure comparable to thermal pressure); 

• I < t^r > I ^ (mean radial speed comparable to the speed of sound); 

• Br,s > \B^^s\ (with Bj. increasing already at z = 0). 

Weakly coupled disks 

• vao C (midplane magnetic field is highly subthermal); 

• \ < Vr > \ C (mean inflow speed is highly subsonic);''^ 

• Br^s < 1^0, s| (with Br taking off only when A increases above ~ 1).^^ 

There are two noteworthy features of weakly coupled disk solutions. First, 
even though the bulk of the disk volume is nearly inert, the disk possesses 
"active" surface layers where A > 1, from which a disk wind can be launched 
in the presence of a large-scale, ordered field (or in which MRI-induced turbu- 
lence can operate; see Gammie 1996). Second, in magnetically threaded disks 
angular momentum is transported vertically even in regions where A is still 
< 1 and i?r ~ 0. This is because a measurable azimuthal field component can 
already exist in these regions {\Bif,/Br\ can be 1 when A <^ 1; see foot- 
note 33) and wc recall that the z(p stress exerted by the field is oc B^^dB^/dz. 
As a result, matter can continue to accrete even in the nominally inert disk 
regions, which would not be possible if only a small-scale, disordered field 
(i.e., MHD turbulence) were responsible for angular momentum transport 
(although in practice turbulent mixing of charges and field might enable ac- 
cretion in these regions also in that case; e.g., Fleming & Stone 2003; Turner 
et al. 2007). This could have implications to the ongoing debate about the 
nature of "dead zones" in protostellar disks (e.g., Pudritz et al. 2007; see 
Chap. V). 

4.5 Stability Considerations 

It was noted in Sect. 4.4 that wind-driving disks should be stable to the most 
rapidly growing linear mode of the MRI. On the other hand, by combining the 

These solutions can thus evade the short-lifetime criticism made against strongly 
coupled disk models, namely, that in the absence of a vigorous mass supply 
(which only happens during the early evolutionary phases of protostars) they 
would empty out on a relatively short time scale (e.g., Shu et al. 2008). Recall 
in this connection that the disk formation models described in Sect. 3.4 yielded 
asymptotic (r — > 0) solutions that correspond to weakly coupled disks. The short- 
lifetime conundrum could be avoided altogether even in a strongly coupled disk 
if angular momentum transport by a wind only dominated in the innermost disk 
regions, which is not inconsistent with existing observational data. 
Using the hydrostatic approximation one can derive a differential equation re- 
lating Br and B^, which, assuming = and a vanishing Hall term, takes the 
form 4^ « —2A at the midplane. 
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wind-launching condition (6) with the condition for the onset of a radial inter- 
change instability (e.g., Spruit et al. 1995), one can show (see Konigl & Wardle 
1996) that, to be unstable to radial interchange, the magnetic term must be 
comparable to the gravitational term in the radial momentum equation. This 
is clearly not the case in the derived disk solutions, which are characterized 
by an inherently subthermal magnetic field (see, e.g., the inequality (2) in 
Eq. (84)). It is thus seen that, when a large-scale magnetic field is responsibe 
for the entire angular momentum removal from an accretion disk through 
a CDW, it automatically lies in a "stability window" in which it is strong 
enough not to be affected by the MRI but not so strong as to be subject to 
radial interchange. 

There are, of course, other potential instabilities to which such disks might 
be susceptible (see Konigl & Wardle 1996 for some examples), but here we 
focus on the question of whether there might be an inherent aspect of the 
wind-related angular momentum transport that could render the disk/wind 
system unstable.^* Based on approximate equilibrium models, Lubow et al. 
(1994b) and Cao & Spruit (2002) suggested that an inherent instability of 
this sort may, in fact, exist (see also Campbell 2009). They attributed this 
instability to the sensitivity of the outflowing mass flux to changes in the field- 
line inclination at the disk surface (^s), according to the following feedback 
loop: 

\vr\ increases ^ tan^g increases 

the wind is loaded by higher — density gas 

(and the removed angular momentum My, I) increase 
\vr\ increases even more.... 

The issue was reexamined by Konigl (2004), who used the disk/wind model of 
Wardle & Konigl (1993) and appealed to the fact that the stability properties 
generally change at a turning point of the equilibrium curve in the solution 
parameter space. Figure 8 depicts such equilibrium curves, labeled by their 
vertically constant Elsasser-number values A = T, in the Blandford & Payne 
(1982) K — X wind parameter space. The lower branches of these curves end on 
the long dashed curve, below which the outfiows remain sub-Alfvenic, whereas 
the upper branches end on the short-dashed curve, to the right of which 
the surface layers of the disk are super-Keplerian. Any particular solution 
is determined by the value of the disk parameter oq = (2/7/3o)^^^, which 
increases along each A = const curve from its minimum value on the "super- 
Keplerian" boundary (see the inequality (1) in Eq. (84)). It is seen that each 



A distinct, but still relevant, question is whether the wind by itself may be un- 
stable. Numerical simulations that treated the disk as providing fixed boundary 
conditions for the outflow have indicated that the wind should not be disrupted, 
even in the presence of nonaxisymmetric perturbations (e.g., Anderson et al. 2006; 
Pudritz et al. 2007). However, these simulations did not account for the existence 
of a feedback between the wind and the disk. 
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curve exhibits a turning point. This behavior can be understood as follows. 

The vertical hydrostatic equilibrium equation implies 6^ ^ « 2/ag (assuming 
:s> B'^s), and since \B^\ evolves with (see footnote 33), it follows 
that \B^^s\/Bz, and hence the magnitude of the angular momentum that the 
outflow must carry away, increases with decreasing ag. Initially, as oq decreases 
from ~ 1, 6r,s = Br^s/Bz increases rapidly, and the corresponding increase in 
the cylindrical radius of the Alfven surface (the effective lever arm for the back 
torque exerted by the outflow on the disk, which scales as X^^"^) increases 
the value of A and leads to a reduction in the ratio of the mass outflow 
to the mass inflow rates (estimated as M^/M^ « 1/[4(A — 1)]). However, 
as ao continues to decrease, the rate of increase of b^.s declines while that 
of \B^^s\/Bz increases, and eventually the mass outflow rate must start to 
increase (with A going down) to keep up with the angular momentum removal 
requirements. The transition between these two modes of enhanced angular 
momentum transport: predominantly by the lengthening of the lever arm 
(on the lower branch) vs. mainly by a higher mass-loss rate (on the upper 
branch) occurs at the turning point of the solution curve. One can determine 
which of the two solution branches is stable and which is not by considering 
the competition between radial advection and diffusion on the magnetic flux 
evolution (best done in the ao — A plane), in analogy with thermal stability 
arguments that analyze the relative magnitude of heating and cooling on each 
side of the thermal equilibrium curve. In this way one flnds that the upper 
branch of the ciirves in Fig. 8 is unstable: this is the branch along which an 
increase in the angular momentum transport is accomplished through a higher 
mass outflow rate, precisely the behavior invoked in the heuristic instability 
argument reproduced above. Since, however, there is an alternative way of 
increasing the angular momentum transport (namely, through a lengthening 
of the effective lever arm), there is also a stable branch of the solution curve. 
Physically, the reason why a perturbation that reduces the field inclination 
to the disk does not necessarily trigger an instability is that an increase in 
br,s also results in greater field-line tension, which tends to oppose the inward 
poloidal field bending. Whether a given solution branch is stable or not is 
determined by the extent to which this stabilizing effect can overcome the 
destabilizing influence of increased angular momentum removal brought about 
by the field-line bending. 

It is possible to argue (see Konigl 2004) that real protostellar systems 
likely correspond to the stable branches of the solution curves. Global nu- 
merical simulations could, however, shed more light on this question, as they 
have already started to do. For example, Casse & Keppens (2002) solved the 
nonideal MHD equations (assuming axisymmetry and a polytropic gas) and 
demonstrated jet launching by resistive disks. These results are significant 
on two main counts: first, they corroborate the basic picture of a magnet- 



This feature of the solution curves is also found in the corresponding plots in the 
ao — A plane, where now k varies along each A = const curve. 
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Fig. 8. Mapping of the wind-driving disk solutions (for several values of the neutral- 
field coupling parameter A) onto the self-similar wind solution space (defined by 
the values of the mass-loading parameter k and angular-momentum parameter A). 
The field-strength parameter oo increases on moving counterclockwise along a given 
curve. 

ically threaded diffusive disk centrifugally driving a wind that removes the 
bulk of its angular momentum, and, second, they indicate that these config- 
urations are stable inasmuch as they reach a quasi-steady state. Subsequent 
simulations have extended this work while confirming its main conclusions 
(e.g., Kuwabara ct al. 2005; Mchani ct al. 2006; Zanni et al. 2007). Three- 
dimensional simulations are the natural next step. 

5 Disk-St£Lr Magnetic Coupling 

We now turn to the best-established example of a large-scale field in proto- 
stars, namely, the stellar dynamo-generated field, and consider the dynamical 
interaction between this field and the surrounding accretion disk. This is a 
rich topic that has implications to the structure of the innermost disk region, 
the stellar rotation and the way mass reaches the stellar surface, protostellar 
jets, and the observational signatures of the region where most of the gravi- 
tational potential energy is liberated. Owing to the nontrivial topology of the 
stellar (and possibly also disk-generated) magnetic field and to the complex- 
ity of the field-mediated interaction between the star and the disk even for 
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simple flux distributions, our understanding is far from complete, with new 
progress being driven largely by improved numerical simulations. We start 
this section by briefly surveying the relevant observations and then describe 
the basic physical processes that underlie this interaction and the results of 
some of the recent numerical work Relevant additional details are provided in 
Chap. I [Protoplanetary Disk Structure and Evolution]. 

5.1 Phenomenology 

We discuss the observational manifestations of two phenomena in which a 
large-scale, ordered magnetic field has been implicated: the relationship be- 
tween protostellar rotation and accretion disks, and magnetic disk truncation 
and resultant field-channeled accretion. The theoretical framework that re- 
lates these two phenomena is outlined in the next subsection. 

Stellar Rotation and Accretion 

It is now well established (based on ~ 1700 measured rotation periods) that 
about 50% of sun-like protostars (M* ^ 0.4 — 1.2 AIq) undergo significant 
rotational braking during their pre-main-sequence (PMS) contraction to the 
Zero- Age Main Sequence (ZAMS), and that the objects whose specific angular 
momenta decrease with time are essentially the ones that arc already slowly 
rotating at an early (~ IMy) age (e.g., Herbst et al. 2007). A clue to the 
braking mechanism is provided by the finding of a clear correlation between 
having a comparatively long period (> 2 d) and exhibiting an accretion-disk 
signature (near-IR or mid-IR excess). This trend is particularly noticeable 
in higher-mass (> 0.25 Mq) protostars, but it is also found in very-low-mass 
stars and brown dwarfs (in which. liowev(;r, the braking efficiency is evidently 
lower). The connection with disks is supported by the fact that the inferred 
maximum stellar braking times (~ 5 — lOMyr) are comparable to the max- 
imum apparent lifetimes of gaseous accretion disks (which are ^ 1 Myr for 
~ 50% of the stars and ^ 5 Myr for almost all stars). These results imply 
that the dominant braking mechanism in protostars is directly tied to active 
disk accretion. In the absence of such accretion (or after the disks disperse), 
PMS stars nearly conserve specific angular momentum and therefore spin up 
as they contract to the ZAMS. As we discuss in the next subsection, stellar 
magentic fields likely play a key role in the braking process. 

Disk Truncation and Field- Channeled Accretion 

As was originally inferred in the case of accreting magnetized neutron stars 
and white dwarfs, a sufficiently strong protostellar magnetic field can be ex- 
pected to truncate the accretion flow, with the truncation radius increasing 
with and decreasing with Ma (see Eqs. (85) and (86) below). In this pic- 
ture, the intercepted matter "climbs" onto the field lines and is magnetically 
channeled to some finite stellar latitude. By the time the matter reaches the 



54 



Arieh Konigl and Raquel Salmeron 



stellar surface it is streaming along the field lines with a near free-fall speed 
and is therefore stopped in an accretion shock. The basic elements of this sce- 
nario are supported by observations of CTTSs (but also of lower-mass brown 
dwarfs and of higher-mass Herbig Ae/Be stars, although so far these have been 
less well-studied). The main findings have been (see Bouvier et al. 2007): 

• The common occurrence of inverse P Cygni profiles, with redshifted ab- 
sorption reaching several hundred km s~^. 

• Observed hydrogen and Na I line profiles that can be adequately mod- 
eled in this picture; the detection in the UV-ncar-IR spectral range of 
predicted statistical correlations between the line fluxes and the inferred 
mass accretion rate. 

• Observed spectral energy distributions of optical and UV excesses that are 
successfully reproduced by accretion shock models. 

• Periodic visual-flux variations due to "hot spots" (interpreted as non- 
axisymmetric accretion shocks on the stellar surface) that are only ob- 
served in actively accreting CTTSs (but not in the weak-line T Tauri 
stars (WTTSs), which lack the signatures of vigorous accretion). How- 
ever, in contrast to "cool spots" (the analogs of sunspots), which last for 

10^ — 10'^ rotations, the hot-spot periodicity only persists for a few ro- 
tation periods, indicating a highly nonsteady configuration. This inference 
is supported by the detection of high-amplitude irregular flux variations. 

• Strong near-IR "veiling" variability that can be interpreted in the context 
of this picture as arising in the interaction region between the disk and an 
inclined stellar magnetosphere, with possible contributions also from the 
accretion shock and from the shock-irradiated zone in the inner disk. 

Zeeman-broadening measurements in a growing number of CTTSs have 
yielded an intensity-averaged mean surface magnetic field strength of ~ 
2.5 kC,^^ with the field inferred to reach ~ 4 — 6kG in some regions. Further- 
more, circular polarization measurements in lines associated predominantly 
with the accretion shock have revealed rotational modulation and demon- 
strated that the field is highly organized (with peak value ~ 2.5 kG) in the 
shock region (covering < 5% of the stellar surface). These results are con- 
sistent with the model predictions. No net polarization has, however, been 
found in photospheric absorption lines, implying a complicated surface field 
topology (with a global dipole component < 0.1 kC). A likely physical picture 
is that there exists a spectrum of magnetic fliix loops of various lengths that 
extend from the stellar surface, with a fraction reaching to ^ 5 — 20 i?* and 
intercepting the disk (see Feigelson et al. 2007). This picture is supported by 
analyses of X-ray flares from young stars (e.g., Favata et al. 2005) as well as 
by direct radio imaging of the large-scale magnetic structures (e.g., Loinard 



A similar value has now been inferred also in a protostar that is at an earlier 
(Class I) evolutionary pliase (Johns-Krull et al. 2009). 
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et al. 2005) and by spectropolarimetric reconstructions of the magentospheric 
topology (e.g., Donati ct al. 2007). 

Numerical studies of magnetospheric accretion indicate that disk mass- 
loaded outflows (which are often predicted to be nonsteady) could be pro- 
duced in the course of the star disk interaction (see Sect. 5.2). There is now 
observational evidence for accretion-induced winds in CTTSs emanating from 
both the inner disk and the star (e.g., Edwards et al. 2006). The stellar com- 
ponent has been inferred to move on radial trajectories and to undergo full 
acceleration (up to a few hundred kms~^) in the stellar vicinity. Although 
this component may well be launched along stellar magnetic field lines, the 
absence of evidence for it in WTTSs suggests that it, too, is triggered by the 
interaction with the disk. 



5.2 Disk-Stcir Coupling Models 

Basic Concepts 

Spherical accretion with Ma — AttR"^ g{R)\V}i{R)\ at the free-fall speed 
vs{R) = (2GM»/i?)^/^ onto a star of mass will be stopped by the mag- 
netic stresses of the stellar magnetic field at a distance where the ram pressure 
of the flow {gv'^) becomes comparable to the magnetic pressure of the stellar 
field (B^/Stt). For simplicity, we assume that the field at that location can be 
approximated by its dipolar component (corresponding to a dipole moment 
yu*). For accretion in the equatorial plane, B{r) = pb^jr^ = B^R^/r^, and one 
obtains the nominal accretion Alfven radius 

4/7 

= ^* , • 2/7 • (85) 

(2GM,)V7Mf/' 

Although the radial ram pressure of the flow is relatively small for disk ac- 
cretion, one can still define in this case a magnetospheric boundary radius 
from the requirement that the total material and magnetic stresses (or energy 
densities) be comparable, i.e., B^/Stt = gv'^ -|- P w gv^. By approximating 
V(j, ~ UK, one obtains an expression similar to the one given by Eq. (85). This 
radius determines the region within which the stellar magnetic field controls 
the flow dynamics and provides a lower bound on the inner disk radius. To es- 
timate the radius where matter leaves the disk, note that the physical process 
of disk truncation requires that the torque exerted by the stellar field that 
penetrates the disk be comparable to the rate at which angular momentum is 
transported inward by the accretion flow, —r'^B^B^ s = Afad(r^J7)/dr. Taking 
|B<^,s| ~ Bz (where we use the approximate maximum value of the azimuthal 
field component at the disk surface; see Uzdensky et al. 2002a and cf. Eq. (38)) 
and setting |d/2/dr| ^ f^K/r yields an estimate of the disk truncation radius 
Td that is of the order of ta, 



rd = kA r A 



(86) 
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Semianalytic models (e.g., Ghosh & Lamb 1979) and numerical simulations 
(e.g., Long at al. 2005) yield kp^ « 0.5. 

If the star is rotating with angular velocity 12, , one can define the corota- 
tion radius by setting OK{r) = f^*, which gives 



(87) 



The interaction of the disk with the stellar magnetic field naturally divides 
into two qualitatively different regimes depending on the ratio r^/rco- 

Td < rco funnel— flow regime 

''d ^ ^co propeller regime . 

In the funnel-flow regime, the disk angular velocity is higher than that of the 

star, and matter can climb onto the stellar magnetic field lines that thread 
the disk and reach the stellar surface. By contrast, in the propeller regime the 
stellar angular velocity exceeds i7K(''d): when disk matter becomes attached 
to the stellar field lines its angular momentum increases above the rotational 
equilibrium value for that radius and it moves outward. Numerical simulations 
have found that, if the disk effective viscosity and magnetic diffusivity are 
relatively high, most of the incoming matter is expelled from the system, 
both in the form of a wide-angle CDW from the inner regions of the disk 
and as a strong, collimated, magnetically dominated outflow along the open 
stellar field lines near the axis (e.g., Romanova et al. 2005; Ustyugova ct al. 
2006).'^^ This mechanism could be relevant to the initial spindown (on a time 
scale < 10*^ yr) of CTTSs. 

Funnel-Flow Regime 

Ghosh & Lamb (1979) proposed that a "disk locked" state in which the torque 
exerted by the field lines that thread the disk (as well as by the material that 
reaches the star) could keep the star rotating in equilibrium (neither spinning 
up nor down). In this picture, the field lines that connect to the disk within the 
corotation radius (as well as the accreted material) tend to spin the star up, 
whereas the field lines that connect at r > r^o have the opposite effect. This 
scenario could potentially explain the relatively low rotation rates observed 
in CTTSs (e.g., Konigl 1991).38 

This picture was, however, challenged on the grounds that the twisting 
of the magnetic field lines that thread the disk at r > Tco will tend to open 



A similar outflow pattern, comprising a conical disk wind and a higher-velocity, 
low-density jet component, was found also in simulations of the funnel-flow regime 
(e.g., Romanova et al. 2009). 

Naively, one would expect a star accreting from a rotationally supported disk to 
rotate near breakup speed; in reality, CTTSs rotate, on average, at about a tenth 
of this speed. 
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them up, thereby reducing the spin-down torque on the star by more than 
an order of magnitude compared to the original calculation (Matt & Pudritz 
2004, 2005). This argument is based on the fact that stellar field lines that 
connect to the disk at r 7^ rco are twisted by the differential rotation between 
their respective footpoints. Initially, the twisting leads to an increase in |i30,s|, 
but then the built-up magnetic stress causes the field lines to rapidly elongate 
in a direction making an angle ~ 60° to the rotation axis (e.g., Lynden-Bell & 
Boily 1994). During this phase |i?0,s| decreases as the field- line twist travels 
out to the apex of the elongating field line, where the field is weakest, a process 
that can be understood in terms of torque balance along the field line (e.g., 
Parker 1979). As the twist angle approaches a certain critical value (^ 4rad 
for a star linked to a Keplerian disk), the expansion accelerates and the mag- 
netic field formally reaches a singular state (a "finite time singularity"), which 
in practice means that it opens up. While the twisting can be countered by 
magnetic diffusivity in the disk, this is only possible if the steady-state surface 
azimuthal field amplitude {2TrrSAf2/ 0^)82,3, where S is the vertically inte- 
grated electrical conductivity (treated as a scalar) and Ai7 is the differential 
rotation rate, does not exceed the maximum value of \B^^s\ in the absence of 
diffusivity (e.g., Uzdensky et al. 2002a). The conductivity in the innermost 
regions of protostcUar disks is probably too high for a steady state to be feasi- 
ble except in the immediate vicinity of Tco (see also Zweibel et al. 2006), and 
in some proposed models (e.g., the X-wind; see Shu et al. 2000) it has, in fact, 
been postulated that the stellar field lines can effectively couple to the disk 
only in that narrow region. 

Numerical simulations have, however, verified that, when the star has a 
strong enough field to disrupt the disk at a distance of a few stellar radii and 
channel accreting matter along field lines, it can maintain an equilibrium disk- 
locking state in which i?* is close to the disk angular velocity at the truncation 
radius (Long et al. 2005). In these simulations it was foimd that, in equilib- 
rium, rco/^d ~ 1.2 — 1.5, very close to the prediction of the Ghosh & Lamb 
(1979) model. Furthermore, as envisioned in the latter model, closed magnetic 
field lines that link the star and the disk exert the dominant stresses in this 
interaction. However, in contrast to the original picture, this linkage does in 
fact occur primarily near Tco- The torque balance is achieved in part through 
field-line stretching (which mimics the connection to material at r > r^-o in 
the original model) and by magnetically driven outflows (which, however, are 
found to remain comparatively weak in the simulations). Field-line opening is 
not a major impediment to this process, in part because opened field lines tend 
to reconnect (especially if the departure from axisymmetry is not large). 



The idea that twisted field lines that open up could subsequently reconnect, 
leading to a repetitive cycle of inflation and reconnection and resulting in the 
star-disk linkage being steady only in a time-averaged sense, was first proposed 
by Aly & Kuijpers (1990) and received support from subsequent investigations 
(e.g., Uzdensky et al. 2002b; Romanova et al. 2002). 
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The question of whether (or under what circumstances) the magnetic 
transfer of angular momentum to the disk is the most efficient way of attaining 
protostellar spin equiUbrium is, however, still unresolved. For example. Matt 
& Pudritz (2008) suggested that stellar winds driven along open magnetic 
field lines could dominate the braking torque on the star. These winds are 
inferred to be powered by the accretion process, but it remains to be deter- 
mined how this could happen in practice and whether the proposed outflows 
are, in fact, related to the stellar winds already identified observationally. It 
is also important to bear in mind that the star-disk magnetic linkage mech- 
anism may, in reality, be more complex than the simplified picture outlined 
above. Some of these expected complications are already being investigated 
with the help of 3D MHD codes, including the effects of misalignment between 
the rotation and magnetic axes, of an off'-centered dipole, and of higher-order 
multipole moments (e.g., Romanova et al. 2003, 2004, 2008; Long ct al. 2007, 
2008; Kulkarni & Romanova 2009).^° Furthermore, the field topology and the 
nature of the interaction could be modified if the disk itself contains a large- 
scale magnetic field (e.g., Hirosc et al. 1997; Miller & Stone 2000; Ferreira et 
al. 2000, 2006; von Rekowski & Brandenburg 2004, 2006). 

Nonsteady accretion 

The magnetic interaction between stars and disks could be variable on a time 
scale as short as i7~^. One possibility, which could naturally give rise to 
observable "hot spots," is for the system to lack axisymmetry — either because 
of a misalignment between the magnetic (e.g., dipole) and rotation axes (e.g., 
Romanova et al. 2004, 2008; Kulkarni & Romanova 2009) or because of the 
intrinsic structure of the magnetic field (e.g., von Rekowski & Brandenburg 
2004, 2006). A longer variablity time scale is implied by the suggestion of 
Goodson & Winglee (1999) that the truncation radius would oscillate on the 
diffusion time scale of the field into the disk (resulting in field-line reconnection 
events and episodic polar ejections) if the diffusivity at the inner edge of the 
disk were relatively low (so \{vB,r/vr) — 1| < !)• A different type of variability 
is implied by the star/disk dynamo model of von Rekowski & Brandenburg 
(2006). In this model, the magnetic field geometry changes at irregular time 
intervals (with magnetic polarity reversals occurring mostly on time scales of 
less than a day), and the star/disk system alternates between magnetically 
connected and disconnected states. The model predicts strong outflows, both 
from the inner disk and from the stellar surface, but typically only a small 
fraction of the disk accretion flow reaches the stellar surface. In fact, there is an 
anti-correlation between the stellar magnetic field strength and the accretion 
rate, and material that reaches the stellar surface comes in at a low velocity 

Given the importance of the question of how matter crosses magnetic field lines 
in this problem, it would be helpful to examine these eff'ects using an explicitly 
resistive numerical code. So far, however, only the axisymmetric version of this 
scenario has been investigated in this way (e.g., Bessolaz et al. 2008). 
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(von Rekowski & Piskunov 2006). It remains to be determined whether these 
two aspects of the model are consistent with the polarization data and the 
evidence for accretion shocks (see Sect. 5.1). 

6 Conclusion 

The discussion in this chapter can be summarized as follows: 

• There is strong observational evidence for a disk wind connection in pro- 
tostars. Large-scale, ordered magnetic fields have been implicated theoret- 
ically as the most likely driving mechanism of the observed winds and jets. 
The ubiquity of the outflows may be related to the fact that centrifugally 
driven winds (CDWs) are a potentially efficient means of transporting 
angular momentum from the disk. 

• Ordered magnetic fields could arise in protostcllar disks on account of (i) 
advection of interstellar field by the accretion flow, (m) dynamo action in 
the disk, and {Hi) interaction with the stellar magnetic field. 

• Semianalytic MHD models have been able to account for the basic struc- 
ture of diffusive disks that drive CDWs from their surfaces as well as for 
the formation of such systems in the collapse of rotating molecular cloud 
cores. Some of these models already incorporate a realistic disk ionization 
and conductivity structure. These studies have established that vertical an- 
gular momentum transport by a CDW or through torsional Alfven waves 
(magnetic braking) could in principle be the main angular momentum 
removal mechanism in protostellar disks and determined the parameter 
regime where wind transport can be expected to dominate radial trans- 
port by MRI-induced turbulence. Further progress is being made by in- 
creasingly more elaborate numerical simulations (involving nonideal MHD 
codes) that have started to examine the global properties, time evolution, 
and dynamical stability of the magnetic disk/wind system. 

• Robust observational evidence also exists for a magnetic interaction be- 
tween CTTS disks and their respective protostars, including strong indica- 
tions of a field-channeled flow onto the stellar surface. This interaction is 
likely to involve mass ejection and is thought to be responsible for the com- 
paratively low rotation rates of CTTSs. Since the magnetic field geometry 
in the interaction region is evidently quite complex and the interaction is 
likely time dependent, numerical simulations are an indispensable tool in 
the study of this problem. 

Future advances in this area will probably arise from a combination of new 
observational findings, the refinement of ciirrent theoretical approaches, and 
the incorporation of additional physics into the models. On the observational 
side, the main challenge is still to demonstrate the existence of CDWs in 
protostars and to determine their spatial cixtent (spread out over most of the 
disk surface or occurring only near its inner edge, and, if the former, is the 



60 Arieh Konigl and Raquel Salmeron 

launching region nearly continuous or is it confined to localized patches?). 
Recent attempts to measure rotation in the outflows could potentially help 
to answer this question, but, as noted in Sect. 2.3, the results obtained so far 
are still inconclusive. 

Regarding the further development of theoretical tools, the greatest im- 
pact would likely be produced by numerical simulations that study vertical 
angular momentum transport by either a CDW or magnetic braking with 
codes that include a realistic conductivity tensor and that have full 3D and 
mesh-refinement capabilities. Such simulations should be able to clarify the 
relative roles of vertical and radial angular momentum transport and the 
possible interplay between them for relevant combinations of the disk model 
parameters (see Sect. 4.4). An interim step might be to solve for the evolu- 
tion of a vertically integrated disk whose properties, at any radial grid zone, 
arc determined from a vertical integration of a simplified version of the radi- 
ally localized, steady-state disk model described in Sect. 4.4. Time-dependent 
models of this type could examine the behavior of a magnetically threaded 
disk after its mass supply diminishes or stops altogether (corresponding to 
the protostellar system evolving into the optically revealed phase) , which has 
previously been studied only in the context of "a- viscosity" models. As was 
already noted above, a state-of-the-art, 3D, nonideal-MHD code is also crucial 
for investigating the stability of disk/wind systems and the various aspects 
of the star-disk field-mediated interaction. One could, however, also benefit 
from further development of the semianalytic models, which could include an 
extension of the ionization/conductivity scheme, a derivation of self-similar 
disk/wind solutions that allow for a radial drift of the poloidal magnetic field 
(see Sect. 4.4), and a calculation of the predicted observational characteristics 
of wind-driving disks. 

The mass fraction and size distribution of dust grains in the disk have 
a strong effect on its ionization and conductivity structure and on the de- 
gree of field-matter coupling (see Sect. 4.4). Existing models incorporate the 
effect of dust in a somewhat ad-hoc manner, by adopting an assumed distribu- 
tion. In reality, the grain distribution is determined by the balance of several 
processes, including grain collisions due to relative velocities that develop as 
a result of Brownian motion, differential vertical-settling speeds, and turbu- 
lence, which can lead to either coagulation or fragmentation. Grains are also 
subject to a coUisional drag force exerted by the gas and arising from the fact 
that the gas is subject to thermal and magnetic forces that do not affect the 
dust. This leads to vertical settling as well as to radial migration, directed 
either inward or outward depending on whether the gas rotation is sub- or 
super-Keplerian, respectively. Furthermore, radial or vertical gas motions can 
affect sufficiently small grains through advection. Yet another effect is evap- 
oration by the ambient radiation field, which could impact dust located at 
sufficiently high elevations and small radii. Some of these effects have already 
been incorporated into generic viscous disk models (e.g., Dominik et al. 2007; 
Brauer et al. 2008) , and one could similarly consider them in the context of a 
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wind-driving disk model. In view of the fact that the latter model is charac- 
terized by a vertical outflow and by comparatively fast radial inflow speeds, 
one can expect to find new types of behavior in this case. In particular, grains 
located near the disk surfaces would either settle to the midplane if they are 
large enough or be uplifted from the disk if they arc sufliciently small, whereas 
intermediate-size grains would first leave the disk and then re-enter at a po- 
tentially much larger radius, from which they could be advected back inward. 
By including dust dynamics, one could examine whether the effect of dust on 
the gas motion (through its infiuence on the field-matter coupling) and the 
effect of gas on the grain motions (through gas-dust collisions) together place 
meaningful constraints on the resulting grain distribution. One could also in- 
vestigate whether the predicted radial transport of intermediate-size grains 
from small to large radii and their possible thermal processing outside the 
disk could be relevant to the accumulating evidence for an outward transport 
of crystalline grains in the protosolar nebula and in other protostellar disks, 
and whether the implied dust distribution in the disk and the wind might have 
distinct observational signatures that could be tested by spectral and imaging 
techniques (see Millan-Gabet et al. 2007, Sect. 4.3, and Chaps. I and VIII). 

Dust particles are thought to be the building blocks of planetesimals and 
their distribution in the disk is thus a key ingredient of planet formation mod- 
els. In fact, the general properties of a protostellar disk are evidently relevant 
to planet formation in light of the growing evidence that the latter is strongly 
influenced by physical processes that occur when the disk is still predomi- 
nantly gaseous. A disk threaded by a large-scale, ordered magnetic field could 
potentially have unique effects on planet growth and migration. One such ef- 
fect is the generation (through magnetic resonances that are the analogs of 
Lindblad resonances) of a global torque that may reduce, or even reverse, the 
secular inward drift (the so-called Type I migration) predicted for low-mass 
planets, which has posed a conundrum for current theories of planet forma- 
tion. As was demonstrated by Terquem (2003) and Fromang et al. (2005), a 
torque of this type could be produced if the disk had a comparatively strong 
(MRI-stable, but still subthermal) azimuthal field (with a nonzero vertical av- 
erage of B^) that fell off sufficiently fast with radius (oc — r^^). A poloidal 
field component could in principle also contribute to this process (Muto et al. 
2008) . Given that a large-scale field with precisely these properties is expected 
in wind-driving protostellar disks (see Sects. 3 and 4), this possibility clearly 
merits an explicit investigation in the context of the disk models considered 
in this chapter. The influence of the vertical channel of angular momentum 
transport and of the overall effect of an ordered, large-scale field on the disk 
structure in such systems (e.g., the reduction of the density scale height by 
magnetic squeezing) may also be worth examining in this connection. 
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